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Abstract: A simple holographic model is presented and analyzed that describes chiral symmetry 
breaking and the physics of the meson sector in QCD. This is a bottom-up model that incorporates 
string theory ingredients like tachyon condensation which is expected to be the main manifesta- 
tion of chiral symmetry breaking in the holographic context. As a model for glue the Kuperstein- 
Sonnenschein background is used. The structure of the flavor vacuum is analyzed in the quenched 
;> approximation. Chiral symmetry breaking is shown at zero temperature. Above the deconfmement 
^ ! transition chiral symmetry is restored. A complete holographic renormalization is performed and the 
^ chiral condensate is calculated for different quark masses both at zero and non-zero temperatures. 
'The ++ ,0 h , 1 ++ , 1 meson trajectories are analyzed and their masses and decay constants are 
computed. The asymptotic trajectories are linear. The model has one phenomenological parame- 
ter beyond those of QCD that affects the 1 ++ , h sectors. Fitting this parameter we obtain very 
good agreement with data. The model improves in several ways the popular hard-wall and soft wall 
bottom-up models. 
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1. Introduction 

The AdS / CFT correspondence |l| has been one of the most fruitful arenas for research in fundamental 
physics in the last decade. Having the possibility of mapping strongly-coupled field theories to weakly 
coupled gravity has set the stage for a large amount of effort devoted to use this idea in order to 
obtain new results on the physics of the strong interactions (see p[, for an introduction). It is clear 
that a precise and controllable string dual of QCD is far from our present understanding. However, 
it is possible to build models that describe interesting strong coupling phenomena which share many 
similarities with real world physics. For the physics of pure glue, models descending from string 
theory, |3], ||, |J have provided important clues towards the confinement of color. Phenomenological 
models for glue inspired and motivated from the AdS / CFT correspondence ranged from very simple 
like AdS/QCD || to more sophisticated versions, namely "Improved Holographic QCD" that capture 
rather accurately the dynamics of glue at both zero and finite temperature ||. 

In this context, the holographic description of chiral symmetry breaking (xSB) has been thor- 
oughly studied. Early examples were studied in the Maldacena-Nunez 0] and Klebanov-Strassler 
backgrounds ||. Later, more QCD-like xSB, in the sense that the operator condensing was bilinear 
of fundamental fields, was described in 0, In the beginning chiral symmetry breaking involved 
abelian chiral symmetry, ||, [5], || [10| . A major breakthrough was the Sakai-Sugimoto model [|l 1 
where the broken symmetry is non-abelian U(Nf)i x U(Nf)n — > U(Nf)v, as opposed to just a U(l). 



However, the model of [11] and generalizations of it have its own limitations. Just as an example, 



one can mention the absence of a tower of excited pions. 
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All the aforementioned models come from well controlled approximations to string theory and, 
accordingly, are top-down approaches. A different possibility is to use string theory just as inspiration 
and define ad hoc holographic models using some QCD features as inputs. This is the bottom-up 
approach. The obvious question is whether the output extracted from such models is larger than their 
input, and we believe that experience has shown that the answer is yes. The benchmark bottom-up 
models for chiral symmetry breaking and meson physics are the hard wall model |T2l |TB| based on 
the Polchinski-Strassler background H] and the soft wall model |T1 . 

In the present work, we analyze in detail several aspects of a bottom-up model presented first in 
It gives a explicit realization of the framework for chiral symmetry breaking first advocated in 
[|l(J: the quarks and antiquarks are introduced by a brane-antibrane pair. A key point of the dynamics 
is the condensation of the lowest lying bifundamental scalar which comes from the strings connecting 
the brane to the antibrane. Around flat space this scalar has negative mass squared and this is 
called the tachyon. This corresponds precisely to a QCD-like chiral symmetry breaking. We need an 
effective action to describe this dynamics and we will resort to the tachyon-DBI action proposed by 
Sen in flT7|| . Once we have decided that we will use this action, we still have to choose an expression 
for the tachyon potential and a holographic geometry in order to use as curved background. For 
these choices, we will constrain ourselves to the (arguably) simplest possibilities. Another interesting 
property of considering brane-antibrane tachyonic actions is that there is a natural Wess-Zumino term, 
which correctly incorporates to the model features like parity and charge conjugation symmetries and 



anomalies; see the discussions in 16 



So far our ingredients are those of a top-down approach. However, it is not possible to stay in 
a limit in which the approximation to string theory is controlled if one wants to reproduce some 
QCD features in this context. For instance, we will need a background with curvature comparable 
to the string scale. We do not regard this point as a negative feature, but just as a sign that our 
approach should be considered of the bottom-up type. This has its own advantages, since for instance 
it seems impossible to get Regge trajectories for excited mesons m? n ~ n from any top-down approach, 
because in those cases the meson mass scale is parametrically smaller than the QCD string tension 
scale, see for instance fll8| . Therefore, the model discussed in this work should be regarded as a 
phenomenological model, partly inspired by top-down consideration and in particular by Sen's action 
[ |17j1 . These top-down inputs will generate some dynamics (compared to other bottom- up approaches) 
which will be crucial in the successful modeling of several QCD features, see the discussions in [|16| 
and also in section 167 



We will constrain ourselves to the abelian case of a single quark flavor but, unlike |)], [TO], this 
is not an essential limitation, since we can make the model non-abelian by piling up branes and 
antibranes, in the spirit of [[□]]. This elaboration, however, is left for future work. 

In section |^, we will discuss the backgrounds (both for confined and deconfined phases) and the 
gravity action of which they are solutions. In section |3], we will study in detail the equation for the 
tachyon modulus r and its bulk vacuum expectation value. In other words, we find the open string 
vacuum and show how it dynamically breaks the chiral symmetry. In section [|, we discuss in detail 
linearized open string excitations around the vacuum, namely the meson physics. A good review of 
this kind of analysis in different holographic frameworks is [[HI. Apart from remarking several general 
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qualitative properties, we end the section with a quantitative phenomenological analysis. Section |] 
provides a brief analysis of the same kind of excitations, but in the deconfmed phase. We conclude 
in section |6] with several discussions; we will convey the pros and cons of the present model and give 
some ideas for future directions. We have relegated various technical comments to eight appendices. 
In particular, in order to facilitate the reading of the text, we review the meaning of the different 
constants and parameters that will appear throughout the paper and for which physical reasons some 
of them are fixed in appendix [A| . 



2. The gravitational background 



As acknowledged in the introduction, the model we will discuss does not come from any controlled 
approximation to string theory. Notwithstanding, we will follow general insights coming from string 
theory and effective actions developed in that framework, especially in the non-critical setting, 
p4| . In this sense, the meson physics (in the quenched approximation) is described by the dynamics 



of a D4-anti D4 system in a fixed closed string background, J 
We take the following gravitational two-derivative action 



S 



d G x^/gl, 



6) 



\ a 



for the background fields: 

1 1 2 

" 2 6T F ( 6 ) 



(2.1) 



with a constant c. We consider the solution discussed in [2T] whose metric is given by: 

R 2 

ds% = -g tt dt 2 + g zz dz 2 + g xx dx\ + g m dr] 2 = — [dx\ z + f A x dz 2 + f A drf] 



with: 



z 



z 



(2.2) 



(2.3) 



A 



This is the AdS soliton, a double Wick rotation of an AdS@ Schwartszchild black hole. The only 
active RR-form we consider is: 

Qc 



(6) 



x 



(2.4) 



for some constant Q c which is proportional to the number of colors and that will not be important 
in the following. The dilaton is constant: 



(2.5) 



QcV 3 

The coordinate rj is compactified and regularity of the metric at z = z\ requires the following 
periodicity condition: 

AlTX 27T 

r] ~ r] + 6rj , 5r] = —z A = — — . (2.6) 



M. 



KK 



The AdS radius is given by: 



R 2 = ^ 



(2.7) 
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The application of this geometry for a phenomeno logical non-critical strings/gauge duality was first 
discussed in ^(| ^TJ. The solution is dual to 1+4 dimensional gauge theory compactified in a circle 
with (susy-breaking) antiperiodic boundary conditions for the fermions. Thus, the low energy theory 
is 1+3 dimensional confining gauge theory coupled to a set of massive Kaluza-Klein fields. 

One can consider the theory at non-zero temperature by compactifying to Euclidean time t E - 
When both circles t E and 77 are compactified, there is a second solution competing with ( |2.2| ): 

d2 

ds% = — [~f T dt 2 + dx\ + f^dz 2 + dri 2 } (2.8) 

z 

while Q, (p|), O still hold. We have introduced: 

/r = 1 - 4 ( 2 - 9 ) 
and is related to the period of the euclidean time and therefore to the temperature as: 

47T 1 

t E ~ t E + 5*^ , <ft E = — z T = - . 2.10 

5 i 



Since when we Euclideanize ( |2.2|) , (|2.8|) both solutions are related by the interchange t E i], z T -H- z\, 
the symmetry makes obvious that there is a deconfining first order phase transition at 

T c =^ = ^ (2.11) 

27T 47T 



For T < T c , the confining solution (|2.2|) is preferred and, conversely, ( |2.8|) dominates for T > T c . Of 
course, this discussion is just a straightforward generalization of 



3. The tachyon vacuum expectation value 

Our main interest will be to study a "tachyon-DBI" action for a single brane-antibrane pair of the 
form advocated in [XT]. In section we will comment about the literature related to effective 
actions including open string tachyon fields and the possible impact of different choices of actions in 
a holographic model of this kind. 

We take the brane-antibrane pair to be at a fixed value of r\ and we will not consider oscillations of 
the transverse scalar, which has no QCD counterpart 1 . The brane and antibrane are at zero distance 
and are therefore overlapping. We have therefore a 5D model for the quarks embedded in a 6D model 
for the glue. The Sen action reads: 

S = - J d A xdzV{\T\) (V- det A L + y/- det A*) (3.1) 



1 A different construction involving D4-anti D4 in this background was considered in 24 . The present scenario is 
more successful in describing different features of QCD. 
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The quantities inside the square roots are defined as: 

A W Miv = 9mn + —F$ N + ra'X {{D M Tf{D N T) + (D N T)*(D M T)) (3.2) 
9v 

where (i) = L,R and the complex tachyon will be denoted T = re td . Indices M,N run over the 5 
world- volume dimensions while we will use /i, v for the Minkowski directions (indices to be contracted 



using rj^ v ). With respect to p7| , we have included two constants gy> A in (|3.2j ), which are related to 



the normalization of the fields to be discussed later. 

The covariant derivative of the tachyon field is defined as: 

D M T = {d M + iA L M - xA^AT (3.3) 

For the tachyon potential we take: 

V = K,e~^ lT2 (3.4) 

where K, is a constant 2 which in principle should be related to the tension of the D4-branes. The 



gaussian is a simple choice that has been discussed in different situations for instance in [fjg, [27], [28| , 
but we warn the reader that it is not at all top-down derived for the present situation and thus 
should be considered as an ingredient of the bottom-up approach. We will comment further in 
section For book-keeping, let us enumerate here the constants that have been introduced up 
to now: R,a',\,gv,JC,fi,z^. In the following, we will impose, on physical grounds, some relations 
among these constants and in appendix [A| we will summarize these arguments. 

We must first find the vacuum of the theory. We should set 6, Al, Ar to zero because of Lorentz 
invariance, but r must have non-trivial dynamics, at least in the confined phase, as will be argued 
below. We thus discuss here the function r(z) that defines the vacuum. The corresponding reduced 
action reads: 

S = -2K J d 4 xdze-^ 2T2 g} t gLy/g zz + 2<Ka'\{d z r) 2 (3.5) 
and the corresponding equation of motion: 

9zz \gtt 9xxJ 2 \g tt g xx g zz J \2iTa'\ J 

We want to study this equation in both the confined and deconfined backgrounds of section |2|. 
For this, we need to explicitly substitute the components of the metric of each background, as given 
in section |2|. We will make these studies separately in the following subsections. Before that, since 
the UV of both solutions is identical (up to 0(z 5 )), the analysis of the UV asymptotics of ( |3.6|) is 
the same. We find that the near-boundary limit z — > limit is given in terms of the two integration 
constants as: 

u 2 

t = c x z + ^-clz 3 log z + c 3 z 3 + 0(z 5 ) (3.7) 
o 



2 We have included the constant dilaton in /C, in order to avoid unnecessary cluttering of formulae. 
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In order to find this expansion, we have imposed that: 



^4 = 3. (3.8) 

This enforces that the scalar bifundamental operator dual to the scalar field (which has mas m 2 = 
— ^ 2 /(27ra'A)) has UV dimension 3 matching the dimension of qq in QCD. This is in agreement with 
the usual AdS/CFT rule A(A — 4) = m 2 i? 2 . It is worth stressing that (|3.8|) should be understood 
as a bottom-up condition on the parameters determining the open string data /i, a', A and not on R, 
since in the quenched approximation one should not think of the flavor branes affecting the closed 
string background. 

The asymptotic expansions for r in the confined and deconfined backgrounds start differing at 
order 0(z 6 ). On the other hand, the IR behaviour for both cases is very different, as will be discussed 
below. 

3.1 The confined phase 

Inserting the metric for the confining background ( [2.2|) into ( p.6|) , we obtain the following equation 
of motion for the order parameter: 

.// 4/X 2 z/a /3 I / 3 /a \_' / " . ,2_/2 



T » _ + ( __ + + ^_ + ^ | T = o (3.9) 



Before going on, notice that equation ( |3.9D depends on two constants z\ and \i. However, such 
dependence can be easily reabsorbed by redefining the field and radial coordinate as z — > z = z/za, 
t — > f — [it. The plots in this section will be performed by taking za = 1, /i 2 = it, but it is automatic 
to find the solution for different values of the constants by rescaling as mentioned above. 

According to the discussion of ||16|| , since the background is confining, we must require the tachyon 



to blow up somewhere. Heuristically, one can think of the diverging tachyon as a brane-antibrane 
recombination; if the tachyon were finite until the bottom of the space one would have an open brane 
(and antibrane). In it was argued that this would lead to bulk flavor anomalies that do not 
match those of QCD 3 . The fact that confinement requires brane recombination (and therefore, chiral 
symmetry breaking) is a Coleman- Witten-like theorem ]30| for the present set-up, and it is analogous 



to a similar discussion of [ 3~I ] for the Sakai-Sugimoto model |TT| . The difference is that the realization 
of chiral symmetry breaking in |3lJ is geometrical while here it is driven by the field r. 



Equation ( |3.9| ) only allows the tachyon to diverge at exactly the end of space (the tip of the cigar) 
z = z\, see appendix [B] for details. 

In the IR, generically the two linearly independent solutions behave as a constant and y/z — z\ 
and they are regular at the tip. There is however a one parameter "boundary" family of solutions 
that (1) depends on a single parameter (2) diverges at the tip. This is the solution we should allow 



3 Anomalies in the hard wall model have been discussed in P9| . In that case, appropriate IR boundary conditions 
have to be imposed on the gauge fields in order to get rid of the IR contribution to the gauge variation of the Chern- 
Simons term. In our case, that contribution is killed due to the divergent tachyon. 
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in the IR. If we call the single parameter C then the acceptable IR solution is: 



T 



^2(z A - z) i(2 ™ 1 C n g n (z) 



n=0 



where 



9n{Z) 



1 + V D n , m ( 1 - — ) 



For the first few constants we have 
Co — C ) C% = - 
6381505 



13 



Co 



247 



3110VC 7 
and the first few functions 



? c 5 



6/x 2 C ' ~* 72/i 4 C 3 
276207997 



26975 



103680/i 10 C 9 



1296/^ 6 C 5 

1402840243831 
33592320/i 12 C n 



(3.10) 



(3.11) 



(3.12) 



9 / z 
g (z) = 1 - — 1 - - 
20 V 



ff3 2 



1 - 



1479 
3380 



1 - 



8481 
4940 



1 



396189 
82004 




(3.13) 



(3.14) 



(3.15) 



(3.16) 



As C increases, the radius of convergence of this series increases. 

The condition that the solution should end up in the one parameter family described above is 
our '"regularity condition". It relates the two UV initial conditions, the source (mass) C\ and the 
vev (chiral condensate) C3. This is a dynamical determination of the condensate as a function of the 
mass by the condition t(z = z\) = 00. This relation will be found numerically. 

In practice, one has to solve numerically the equation of motion ( p.9[) arranging the asymptotics 
to be ( |3.7| ) in the UV and ( p.!0[ ) in the IR. One can implement a standard shooting routine whose 
inputs are c% and some UV and IR cutoffs, where the numerical solution is required to match the 
mentioned asymptotics. The value of C3 leading to ( p. 10 ) is the limiting point between a behavior of 
diverging derivative of r and a behaviour where r remains finite everywhere, see figure ||. 

In fact, for fixed c\ there are two values of C3 for which r diverges at z\, since r can diverge to 
+00, for a particular c 3 > 0; or to —00, for a particular c 3 < 0, (we are assuming by convention that 
C\ > 0). However, the C3 < solution is unstable and should be discarded. This can be understood 
by comparing the free energy of both solutions or, alternatively, by realizing that there is a tachyonic 
mode in the pseudoscalar sector. In the massless quark case c\ = 0, both solutions are related 
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by r — > —r and are physically equivalent. They are just related by a rotation in the direction of 
the Goldstone pion, which is exactly massless. This behaviour is completely analogous to the one 
described in \W\. For illustrative purposes, we plot in figure |1] the result of numerically integrating 



( fT9|) and the behaviour of t(z) for different values of c 3 . 
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Figure 1: All the graphs are plotted using za = 1, ^ 2 = vr and c\ = 0.05. The vertical line at z = za = 1 
represents the IR end of space (tip of the cigar). On the left, the solid black line represents a solution with 
C3 Rj 0.3579 for which r diverges at z\. The red dashed line has a too large C3 - in particular, it corresponds 
to C3 = 1 - such that there is a singularity at z = z s where d z r diverges while r stays finite (a behaviour 
of the type r = k% — k,2y/ z s — z, that is unacceptable since the solution stops at z = z s where the energy 
density of the flavor branes diverges). The red dotted line corresponds to C3 = 0.1; this kind of solution 
ought to be discarded because the tachyon stays finite everywhere. The plot in the right is done with the 
same conventions but with negative values of C3 = —0.1, —0.3893, —1. For C3 « —0.3893 there is a solution 
of the differential equation such that r diverges to —00. As explained in the text, this solution is unstable. 
Thus, the physical solution for this particular value of c\ is uniquely determined to be the solid line of the 
graph on the left. 



In figure [| we plot the values of c 3 and C obtained dynamically, as a function of c\. 
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Figure 2: The values of C3 and C determined numerically as a function of c\. In the first plot we portray c 3 
in terms of c\, for c\ < 1. In the second plot we show C3 for a larger range of c\. The third plot depicts the 
constant C entering the IR expansion as a function of c%. Again, we have used z\ = 1, fi 2 = it for the plots. 
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3.2 The deconfined phase 

Inserting the metric Q2.8|) in (|3 



, we obtain the following equation for r in the deconfined phase: 



T" + 



-T 



4 

z 



z 2 h 



2 12 
fl T 







(3.17) 



The IR behaviour of this equation is quite different from the one of ( |3.9| ). First of all r is not allowed 
to diverge at any point. The difference with respect to the confining case is that since there is a 
horizon, one can allow the branes not to recombine as long as they end on the horizon. Then, they 
will not generate any anomaly. 

Still, one has to discard solutions for which r' diverges at some z < zt (with r remaining finite). 
Those solutions yield infinite energy density and are physically inconsistent, just as in the confining 
case. It turns out that this condition uniquely selects a value for C3, for which r reaches the horizon 
at z = zt taking there a finite value, say t\ z=Zt = ct, as shown in the first plot of figure [| 

Thus, we have a one parameter family of physical solutions (which again fix C3 in terms of C\ 
as expected on physical grounds). In a similar fashion to the confined case, by redefining the field 
and radial coordinate as z — > z = z/zt, t — > f = fir, the dependence of the equation on these 
two parameters can be reabsorbed. Near the IR, these solutions read, in terms of the parameter 
t(z t ) = c T : 

3c T/ x 9c T , a , _ 2 2 2 



-(zt-z) 



+ /j, 2 c t )(z t - zy + 



(3.18) 



5z T 200z T 

Once Ci is fixed, C3 and ct are dynamically determined by this IR condition, and their values 
can be found numerically; using a standard shooting technique. Notice that for c\ = 0, the solution 
is simply r = and chiral symmetry is unbroken. We display some plots with numerical results in 
figure ||. 
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Figure 3: Plots corresponding to the deconfined phase. All the graphs are plotted using zt = 1, /u 2 = n. 
For the first plot we have taken c\ = 0.05. The solid line displays the physical solution C3 = —0.0143 
whereas the dashed lines (03 = —0.5 and C3 = 0.5) are unphysical and end with a behaviour of the type 
r = k\ — k2y/z s — z. The second and third plots give the values of C3 and ct determined numerically by 
demanding the correct IR behaviour of the solution, as a function of c\. 



3.3 Holographic renormalization, the quark mass and the quark condensate 

On general AdS/CFT grounds, we expect the integration constants ci, C3 of the UV expansion (|3.7|) 
to be related to the source and vacuum expectation value of the boundary operator associated to the 



-10- 



bulk field r, which is the scalar quark bilinear qq. Namely c\ should be, essentially, the quark mass 
and C3 the quark condensate. In this section, we will make this connection precise. 



As has been pointed out many times -see for instance ||32|| , ||33|| -, in QCD the quark mass runs 
all the way to zero in the far UV, a fact that cannot be matched in a holographic model with AdS 
asymptotics (such that the m q we will define is the UV value, which does not run further). If we 
want to make a phenomenological analysis, the most natural option is to identify the m q of the model 
with the QCD quark mass measured at a scale around 1 or 2 GeV. It is conceivable that this feature 
can be ameliorated by using the tachyon action in a holographic setup which incorporates asymptotic 
freedom, as in "Improved holographic QCD" ]7). 

The quark condensate is defined as: 

(,-,> = (3,9) 

where in order to find S ren we have to follow the procedure of holographic renormalization, see for 
instance [55| . The first step is to regularize the action by placing a UV cut-off at z — e, namely 
S reg = j A £, where we have defined, from (|3.5|): 



C = -2/C e-s" T ^«Vfc + WW (3.20) 

Since we are just concerned with the variation of S reg with respect to m q , we compute the functional 
derivative with respect to r: 



dr dr' J J e dz \ dr' J 



and therefore 



SSna =l U^U^W rUsr%\ (3.21) 



(3.22) 



5t dr' 

We are interested in = j^^r 2 -- In order to compute Jjp, one should take into account that C3 
is a non-trivial function of c\. We find by explicit computation, using the UV expansion 

5S reg _ 5 2 1 2c\ 2 3 0, n 1 3, ,2 , 2 



^ OV (3^ + gC^ 2 loge + 2c 3 - -c^ 2 + - Cl a ci c 3 J (3.23) 

where we have disregarded terms that vanish as e — > 0. We now have to write the appropriate covariant 
counterterms that should be added to S reg in order to define the subtracted action S SU b = S reg + S ct : 

S ct = -ICR [ d 4 xV=l (-\ + ^-r 2 + ^r 4 loge + ^ar 4 ] (3.24) 
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where 7 corresponds to the induced metric at z = e, namely a/— 7 = R A e~ A . We have introduced a 
constant a which captures the scheme dependence of the condensate and reflects an analogous scheme 
dependence in field theory. It will be further discussed in appendix |CJ The renormalized action is 
just S ren = lim^o Ssub- It is now straightforward to find: 

= -(2WO? 3 A) (-Ac 3 + c?/i 2 (l + a)) (3.25) 
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Notice that the term with cid Cl cs in (|3.23|) drops out because there is one with the opposite sign 



in that cancels it. We now want to evaluate the quark condensate (|3.19| ). The quark mass is 



proportional to ci, and we take it to be 

m q = (3 Cl (3.26) 

where (3 is a constant. 

The arbitrariness of this multiplicative constant related to the normalization of the fields has 
been stressed (in analogous situations) in |33|, pTfl . We finally obtain 



(qq) = -(2ira'lCR 3 \) 



/r(l + a) 



(3.27) 



3.4 The jump of the condensate at the phase transition 

The first term of the expression for the quark condensate Q3.27| ) depends on the quantity c 3 that 
is determined dynamically via the numerical integration. The second term depends on the quark 
mass and a scheme dependent constant a. We now compute an observable which is independent of 
this second term by finding the jump of the quark condensate when the theory is heated such that 
it undergoes the deconfmement phase transition. Concretely, we take a fixed mass (fixed ci) and 
compare C3 for a confined theory and deconflned theories, such that z\ = zt, namely at the phase 
transition point. We have that A(qq) = (qq) C onf — (qq)deconf = — 4:^(2ira' )CR 3 X) AC3. In figure f| we 
plot AC3, which in practice is nothing else that the difference between the first plot in figure |2| and 
the second plot of figure ^. It turns out to be a monotonously decreasing function, at least in the 
range of c\ which we have been able to study numerically. We plot the result in figure [2[ 



^3 

0.35: \ 

0.30 

0.25 

0.20 

0.15 

0.10 

0.05 



dAc 3 
dci 



C l -0.2 




C] 



Figure 4: The finite jump of the quark condensate and its derivative with respect to c\ when the connnement- 
deconfinement transition takes place. The values z\ = zt = 1 and [i 2 = tt have been used in the plot. 



Let us now discuss how the quark condensate changes when tuning the temperature, while keeping 
fixed the quark mass and the QCD scale. We will plot the quantity: 



R 



(qq)o) 



(3.28) 



where (qq)r is the condensate evaluated at temperature T. We have included the power of T c in the 
denominator in order to make the quotient dimensionless. Let us start by computing the explicit 
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value of (qq)o from Q3.27 ). We will consider small quark masses (compared to the QCD scale or 



Mkk) so we can neglect the last term of ( |3.27|) and use the value C3 ~ 0.37z A computed for small c\ 



in the first plot 4 of figure 0. Inserting the value of z\ in terms of T c (|2.11|) and advancing the value 
of the normalization constant that will be found in ( [4.34 ), we have (qq)o ~ —0.3N c /3Tj?. Turning on 
the temperature but staying below the phase transition, the functions of the metric do not depend 
on the temperature and therefore (qq)R = for T < T c . This is just a consequence of large- N volume 
independence. In order to compute the result in the deconfined phase, we would like to use the values 
of C3 as a function of c\ plotted in figure || From the figure, one can fit, for small c\Zt the value 
of c 3 to be c 3 Zy ~ — 0.286ciZ;r- Using this expression, together with ( |3.26|) , ( p. 10] ), (|4.34|) , we have 
(qq)r ~ 0.09N c m q T 2 . Finally, we reach the result: 



{qq)i 



iV c ^|(0.3/3T c 3 + 0.09m g T 2 ) 



(T>T C 



(3.29) 



We illustrate the behaviour of (qq)n in figure |5[ Notice that, since we are considering light quarks, 
the constant term is the largest until T 3> T c . 



<qq> R 

0.025 



0.020 

0.015 
0.010 
0.005 
0.000 



m T/T c 



Figure 5: Behaviour of {qq)n as a function of the temperature. We have taken N c = 3, /3 = 1, m q /T c = 1/40 
for the plot. 



This is in rough qualitative agreement with lattice results, see figure 4 of 
jump at the phase transition is sharp due to the large iV limit. 



In our case, the 



4. Meson excitations: the confined phase 

Up to now, we have discussed the vacuum (saddle point) of the model. We will now discuss in turn 
the different excitation modes, by expanding the action ( |3.1[ ) up to quadratic order in all the fields. 
In this section, we will only refer to the confined phase and therefore (r)(z) is computed as in section 
31. 

We define the vector and axial vector fields as: 

vm = , A M = (.4.1) 

4 We remind the reader that the plots were done by fixing z\ = 1 and Zt = 1 respectively. The values for generic 
z\, zt are obtained just by rescaling C\ — {c^z^ijz^ 1 and C3 = (c3\ ZA= i)z^ 3 , and similarly in the deconfined phase, 
substituting z^ by zj~. 
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The notation for the associated field strengths will be Vmn, Amn- We use a gauge A z = V z = 0. We 
split the relevant fields as: 

V^x^z) = ^ v (z)V,(x^ , 

A^, z) = z) + A\{x^ z) = ip A (z) ApixT) - <p(z) 9,(P(/)) , 

6(x»,z) = 2${z)V{x») , 

T {x*,z) = {t)(z) + s(x^z) = (r)(z) + ip s (z)S(x^) . (4.2) 

where V M and A^ are transverse vectors d^V p = d^A p = 0. A few comments are in order: we have 
used the residual gauge freedom to make V p transverse. We have anticipated the behaviour of the 
equations of motion in order to write down the terms containing V(x), associated to the pseudoscalars. 
The symbol (t)(z) represents the tachyon vev in the bulk, as discussed in section |[ 

The different bulk fields are dual to the field theory quark bilinears due to the boundary cou- 
plings 5 : f d^xV^Jfy, f d^xA^Jji, f d^xVJ'p, f d A xSJs, where the J7"'s are the different bilinear quark 
currents: J$ = qj^q, J% = q^ 5 q, J P = q^ 5 q, J s = qq - (qq). 

We also define the useful quantity: 

9 ZZ = 9zz + 2na'X{d z (r)) 2 (4.3) 

In the rest of this section, we will discuss the explicit prescriptions to compute the masses and 
decay constants for the different mesonic modes. In particular, the decay constants will require 
computing two-point correlators for which one has to holographically renormalize. We give here the 
complete set of counterterms which make the on-shell action finite. 

Set = -ICR I d A x^\ -\ + ^-r 2 + ^r 4 loge + ^ar 4 + 



2 3 18 ° 12 
+ { ^P^lr P l u5 (V, u V p5 + A, u A p5 )(\oge + ~) + ^fr v (D,Ty(D u T)(loge+ \) | (4.4) 

This expression completes (|3.24j) by including all the active fields we are considering. The terms of 1 /2 
inside the brackets of the second line are finite contact terms that have been chosen for convenience. 
We now discuss in turn each of the modes. 

4.1 Vector mesons 

The quadratic action corresponding to the vector mesons that comes from expanding fl3.1|) reads: 



r\2 



S v = - ' K, I d*xdze-z» 
9v 



1 ,,2^2 



1..A 
2 



ghV^V^ + g xx ~gJd z V^d z V^ 



(4.5) 



where we have constrained ourselves to the confining phase in which g tt = g xx . Here and in the 
following, it should be understood that the //, v indices are contracted using the flat Minkowski 



3 The various discrete symmetries and their realization are detailed in [fT™ 
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metric, since we have explicitly written the factor of g xx = R 2 /z 2 . The equation of motion can be 
easily derived: 

- 1 - dz (e-^ 2r2 g xx gJd^ v (z)) - q 2 ^ v (z) = (4.6) 
e—^gh V ' 

where we have gone to Fourier space and defined the 4d-momentum such that for the eigenmodes it 
corresponds to the mass eigenvalues q 2 = —m v . The above equation explicitly depends on only two 
parameters za and /i 2 . It is easy to check that z\ just gives an overall scale to m v (and, in fact, to all 
dimensionful quantities that will appear later) and fi 2 only enters through the combination f 2 = fi 2 r 2 . 
This was the same combination in the tachyon equation (see the comment below Q3.9D ), and in fact 
one can fix the value of fi 2 without any loss of generality. From now on, we will set fi 2 = it, z\ = 1 
in all the plots, although we will keep the parameters explicit in the equations. 

Finally, notice that ( |4.6| ) depends implicitly on c\ (the quark mass) through the bulk vacuum 
expectation value of r. In short, the vector spectrum given by the model depends just on a multi- 
plicative constant z\ and the parameter c\, namely the quark mass. All the other parameters that 
we have defined drop out from this computation. 



4.1.1 Schrodinger formalism and the mass spectrum 



In order to gain some insight in the problem, let us transform equation ( f4.6|) to a Schrodinger problem, 
following appendix [D[ We immediately read C(z) = M(z) = and: 



A(z) = e"*" T g xx ~gj , B(z) = e"^ T ~gl , (4.7) 
such that the Schrodinger radial variable is defined by: 

„ = /\fe = /\&,. (4.8) 



Notice that u G [0, oo). It is now a straightforward exercise to obtain the Schrodinger-like potential 
( p.7| ), for a given c\. One has to compute numerically r{z) as in section [D], then evaluate ( |D.7[ ) and 
finally implement the variable change ( [4.81 ). Some examples are plotted in figure []. 

We observe that the potentials move up as we increase c%. This is of course expected on general 
grounds, since meson masses should grow with increasing quark masses, but this feature is missing 
from the hard wall or soft wall models. In ||15||, we made a phenomenological fit including the strange- 



strange meson masses in the analysis, finding good agreement with experimental data. 

One can check that the leading contribution to V(u) near u = is of the form V(u) = \u~ 2 + . . . . 
This just comes from the UV AdS asymptotics. Let us now find the leading IR contribution to 
V{u). For large u (namely near z = Z\), we have that g xx ~ R 2 /z\ and g zz ~ 2ira'\(d z T) 2 . Using the 
expressions in appendix 0, a little algebra shows that for large u, we have ^ ~ -^z\d z r and therefore 

1 q 2 

u ~ -^z A r. The function 5 behaves as S w ( -j^J e 4z a what finally leads to V{u) = -^u 2 + 0(u). 
Therefore, V(u) grows quadratically at large u, a fact that leads to standard Regge trajectories for 
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Figure 6: The Schrodinger potential associated to the vector excitation for different values of the quark 
mass. From bottom to top c\ = 0,0.5,1,2,3,4. For illustrative purposes, on the right we plot the second 
normalizable "wavefunction" for c\ = 1. 



large excitation number n |I4"|1 . Asymptotically, the slope of these trajectories is Hindoo = p-, 



as can be found by evaluating 

By using standard numerical shooting techniques, we have computed the mass spectrum. In 
particular, we have computed the first seven states, changing the quark mass parameter in the range 
< ci < 5. We plot some results in figure |7|. It turns out that for small c\ the growth of meson mass 
on the quark mass is linear. This is just what one expects from a Taylor expansion if we consider the 
meson masses as function of the quark masses. This was a result already found in j0|. We have: 



z A m$ » 1.45 + 0.718ci , z A m ( y } « 2.64 + 0.594ci , z A 7n$ « 3.45 + 0.581ci , (ci < 1) 
z A m$ » 4.13 + 0.578c! , z A mf « 4.72 + 0.577ci , z A mf « 5.25 + 0.576c x . (4.9) 



At around ci > 1, the graphs start departing from the straight line, as can be seen on the second 
plot in figure ^. 



4.1.2 Current-current correlator and normalization of the action 



We have discussed the vector spectrum, but we are also interested in the decay constants of each state. 
In order to compute them, we have to fix the multiplicative constant associated to the normalization 
of the action associated to the vector modes. We will follow the reasoning of [12, 113 and match the 
correlator HV to the quark bubble perturbative computation at large Euclidean momentum. In fact, 
all the discussion of this subsection is completely parallel to [|12|, 13fl , since it only depends on the 
asymptotically AdS structure. We however repeat the argument in the present notation for the sake 
of clarity. 

The current-current correlator is defined as: 



d A xe i(lx {j,{x)JM) = (v,»q 2 - gVliW) 



(4.10) 
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Figure 7: Vector meson masses as a function of c\ (proportional to the quark mass). On the left, we plot 
the fitted straight lines together with several points computed numerically, for the seven lowest-lying vector 
modes, in the range c\ < 1. On the right, we go to larger values of c\ and, for clarity, only plot the three 
ligthest states. The dashed lines correspond to the linear fits valid for small c\ whereas the solid line are the 
actual values found numerically. 



As usual, we compute it holographically from the on-shell action. Integrating by parts in ( [4.5| ) and 
adding the counterterm from (|4.4j) , we find: 



9 4 v J (2tt) 
JCR 5 (2na' 2 ) f d 4 q 



Sv= y : \ 7 ^e-^g xx g^V,(q,z)d z V^-q,z) 



+ 

, „ ,^r /( ( V .r)r"(-V.r)(l(,gr+^ (4.11) 
9v J l Z7r 

where V M (g, z) = ip v {q, z)Vq {q) , and ip v (q,z) is the solution to ( f4.6|) subject to V(q,e) = 1 and with 
normalizable behaviour in the IR. At small z, the solution for ip v (q, z) can be expanded in terms of 
two integration constants as: 

il) v = bx{q)+ ( h^ + lh (q)q 2 log z J z 2 + . . . (4.12) 



Substituting this expression into Q4.11]) and taking two derivatives with respect to Vo(q), we find 
that: 

n v(9 >) = ^^Ol^i (4 . 13) 

9v 1 

where we have set b\(q) = 1 consistent with the two-point function prescription and the non-trivial 
g 2 -dependence comes through 62(1?), which has to be found by integrating numerically and demanding 
the physical IR behaviour. 

Before entering into numerical integration, we are interested in computing the limiting behaviour 
for Ily for large q 2 . In order to do this, we consider again the equation written in Schrodinger form 
and notice that, for small z: 

u ~ z , a(u) ~ uTiil^ (u) (4-14) 
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The leading large q behaviour is not affected by the details of the Schrodinger potential, so we may 
just approximate it by an expression that interpolates between its UV and IR behaviours, as discussed 
in subsection |4.1.1| . Namely, we can just write: 

-3>+|-^ + cV)« + g 2 a = 0, (4.15) 



Av? 

where we should take c 2 = -rr. However, we will see that the value of c 2 does not matter for the 
normalization we want to make. ( J4.15| ) is nothing else than the soft wall model of The general 
solution of Q4.15D is: 



,2 



a(u) = kx—^U^-, 0, cu 2 ) + k 2 —^L~} % (cu 2 ) (4.16) 

\J u 4c \J u — 



4c 



where U stands for the confluent hypergeometric function and L for a generalized Laguerre polyno- 
mial. IR normalisability requires k 2 = 0. We now substitute in (f4.14| ) and fix ki by demanding that 
\im z ^ il) v (q,z) = \im u ^ u2a(u) = 1. 



* v ^ = i r (i)^ u ^ ^ (417) 



We can now expand this expression for small u ~ z and compare to ft4.12|) , in order to read b 2 and, 
accordingly U v from ( [4.13|) . The leading pieces at large q 2 for b 2 read: 



lim ^| = Ilogg 2 -i(l + log4-2 7 )-^ + ... , ( J. IS) 

g 2 -s>oo q z 4 4 3g 4 



where 7 is Euler's constant. Therefore, the leading piece which we can compare to the quark bubble 
via fl4TT3D is: 

(q 2 ) = K ~ 1 —^ logg 2 4.19 

9v 



By matching this expression to the perturbative result, we find 6 : 

(2W) 2 /CR _ JV^ 
7v "T^ 2 



(4.20) 



One may wonder how good the results obtained from the simple Schrodinger problem we have dis- 
cussed (^4.15|) are as an approximation to the full problem ( [4.6|) . In figure j|, we compare ( |4. 1 8| ) to the 
value of b 2 (q 2 )/q 2 computed numerically. 



6 Notice that we are dealing with abelian flavor symmetry. There is a factor of i difference with respect to jl2[ since 
in that paper they deal with a non-abelian case and define Tr(t a t b ) = \8 ab . This also makes different the definition of 
the decay constants, for instance the defined in [O is the we will use divided by \/2. 
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-0.3 



Figure 8: We plot the value of proportional to the vector current-current correlator. The plot is a 

comparison of the approximation found in the text from a simplified Schrodinger problem, Eq.( 4.1S| ) (red 
dashed line) to the actual numerical result (solid line). The numerical plot was made by taking c\ = 0.1, 



TT. 



4.1.3 Decay constants 

We are now interested in determining the decay constants of our mesonic states. We start by writing 
the current- current correlator as a sum rule. 



(q 2 + — it 



(4.21) 



The idea is to derive the form of the sum rule holographically. In appendix [E] we give a general 
description on how to write holographically a two-point correlator as an infinite sum. Then we use 
properties of the normalizable modes in order to determine the values of F n . The argument follows 
[O, [L3j so we directly quote the result in the present notation 7 : 



2„/,(") 



N c R d z tp 



6n 2 m 2 \ dz 2 



(4.22) 



where i/j- 



(n) 
V ' 



n 



1,2, 



00 are the solutions of 



normalized as: 



B(z)(ipy^) 2 dz — 1 



(4.23) 



with B(z) given in (|4.7|) . Again, we can compute numerically the values of the decay constants given 
by the model. We have plotted them in figure ^|. One can see that the dependence on excitation 
number is rather mild for small quark masses and for a large number of modes. 



7 Notice our definition of F„ is different from [O 
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\n=l 

1 ' ' 2 ' ' 3 ' ' 4 ' : x5 Cl 

Figure 9: The decay constant, in units of z^ 1 for the four lowest-lying, the seventh and the twelve-th vector 
mode (from bottom to top), as a function of c\. The numerical plot was made by taking // 2 = ir and N c = 3. 



4.1.4 Regge trajectories for vector mesons and linear confinement 

Typical holographic models lead to a behaviour of the masses with the excitation number as oc n 2 , 
for large n ||39|| . However, experiment and semiclassical quantization of a hadronic string (assuming 
linear confinement) suggest that vp? n oc n in QCD. Circumventing this problem was the motivation 
for developing the soft- wall model . As pointed out above and also in [HJ, a model including 



an open string tachyon with action ( |3.1| ) and gaussian tachyon potential, naturally implements this 
behaviour. In figure [10], we plot the results of some numerical computations which display this feature. 
We remind the reader that, as we saw in section STTTTj , for vector mesons Hindoo m^ +1 — m 2 = Q/z\. 
This seems to be born out by the figure. 



m„ 
15 



6.2 
6.0 
5.8 
5.6 
5.4 
5.2 
5.0 
4.8 



Figure 10: Results corresponding to the forty lightest vector states with c\ = 0.05 and c\ = 1.5. On the 
right, the horizontal line signals the asymptotic value 6 of the Regge trajectory, the lower line corresponds 
to ci = 0.05 and the upper line to c\ = 1.5. Masses are given in units of z~^ 1 . 
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4.2 Axial-vector mesons 

The quadratic action corresponding to the axial vector mesons that comes from expanding ( |3.1| ) picks 
an extra term with respect to ([4.5|), coming from the covariant derivative of the tachyon: 



(2W)V f \\~ g iA^ + gxx fJd z A,d z A» + ^ 

2 S{27ia' 



S A = - v ~"7 y K, I (fxdze-^ 2 
9v 



—g zz A ilv A ilv + g xx g zz 2 dzA^dzA^ + nfn __ V M ^ 2 t 2 g^gizA^A^ 

(4.24) 



The equation of motion can be derived to be: 



-d z (e-^ 2 g xx g;/d z ^ A (z)) - k^^ A (z) - q 2 ^ A (z) = (4.25) 

where we have introduced a new constant k as the combination: 

k AR ^v (426) 

We observe that r only enters through the combination fi r so \x is immaterial since it can be rescaled 
away. On the other hand, the constant k, that did not enter the parity even sector does affect the 
physics. In fact, by comparing ( (4.6| ) to ( |4.25| ), one can see that the difference between the equation 
for the vectors and the one for the axials in controlled by k. Therefore, it is natural to guess that k 
somehow enhances or suppresses the effects of chiral symmetry breaking on the P-odd spectra. In the 
following, we will see how the physics depends on this parameter. The model of our previous work 
HJ was more constrained since, in terms of the present notation, k was fixed to ■±§. 



4.2.1 Schrodinger formalism and the mass spectrum 

The functions for converting to a Schrodinger problem A(z), B(z) are as before ( (4.7| ). On top of 
that, we have here a non-trivial M(z) given by M(z) = B(z) kfi 2 T 2 /z 2 . It is easy to check that the 
leading piece in the UV of the Schrodinger potential is 3/4u 2 as for the vectors. However, the leading 
IR behaviour is modified due to the term proportional to k to Vjr(u) = (l + 4^) u 2 . From this 
observation, one can immediately realize that the model gives different Regge slopes for vectors and 
axials and that the leading behaviour of ITa^ 2 ) at large Euclidean momentum coincides, consistently, 
with the vector one ( |4.19| ). We will later comment further on these issues. 

The qualitative appearance of the Schrodinger potentials for the axial excitation is similar to 
the ones for the vectors. But the value of the potentials in the axial case is always higher due to 
the terms coming from M(z). Thus, for equal excitation number and quark mass, the axial mode is 
always heavier than the vector mode (with the difference controlled by k). For small values of c\, the 
dependence of the meson masses on the quark masses is linear 

z K m { l ] sa 2.05 + 1.46ci , z A m% ] » 3.47 + 1.24ci , z A mf « 4.54 + 1.17a , (ci < 1) 
z A mf ss 5.44 + 1.13ci, z A mf w 6.23 + l.llci , z A m^ w 6.95 + l.lOci . (4.27) 

For this calculation, we used k = ^| as it is found by the fit of the parameters in section (|4.5|), whereas 



in 



15 1, we used k = For larger ci, the plots of the meson mass dependence on the quark mass for 



the axial excitation look similar to the vector case, see the plot on the right of figure [7]. 
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4.2.2 Current-current correlator and the pion decay constant 

By explicit computation it is easy to check that the UV expansion of the solution to (|4.25|) is given 
in terms of the two integration constants as: 

^ A = h+ f^ + h^q 2 + nkcl) log zj z 2 + ... (4.28) 



For the case of the axial vector, the corresponding sum rule generalising ( |4.21| ) reads: 



2 



Mf) = 4 + E (2 M K 2 ■ \ ( 4 - 29 ) 

Q yT + m n ~~ le ) 

where of course now the n run over the axial resonances. The F n here are computed in essentially 
the same way as for the vector case, namely using Q4.22 ) with a normalization condition ( |4.23| ). 



Now, we are also be able to compute the value of f^. We do this by directly computing the 
2-point function at zero-momentum, namely: 



fl = -^- 2 h\^ ( 4 - 3 °) 



where we have used the expansion (|4.28|) , which up to that order, is also valid for the axial case. The 
value of 62 to be inserted in (|4.30| ) is found numerically by solving fl4.25| ) with q 2 = 0, with initial 



condition ip A \ z=e = 1 and demanding IR normalizability. 





Figure 11: The pion decay constant and its derivative as a function of c\ - the quark mass. The different 
lines correspond to different values of k. From bottom to top (on the right plot, from bottom to top in the 
vertical axis) k = The pion decay constant comes in units of z~j^ ■ 



From the figure, we observe that the decay constant grows with the quark mass for small quark 
masses and then starts decreasing. We can also observe that increasing the parameter k, increases 
the value of the pion decay constant. This is in agreement with the intuitive notion given above that 
k somehow controls the amount of chiral symmetry breaking. 
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4.3 Scalar mesons 

We now deal with the scalar excitation. The quadratic action reads: 



S = -2W/CA / (fxdze—^ T 



_ 3 _ 1 

9 2 xx 9zzgzz 2 (d z s(x, z)) 2 - 2fx 2 g 2 xx g zz 2 r(z)r(zy 's(x, z)d z s(x, z) 

+ ^Tl^ 2 < z f ~ l )9l x ~gis{x )Z ) 2 + g xx g zz fJ{d^s{x ) z)) 2 \ (4.31) 

From ( |3.1| ), it can be seen that there is also a linear term in the bulk action, but can be easily shown 
to be a total derivative. 

We can read the functions that are used to rewrite this problem is Schrodinger form, as defined 
in appendix [D]: 

A I \ --U 2 T 2 2 9zZ n , \ 1 U 2 T 2 9zZ 

A(z) = e a" g z xx ^ , B{z) = e g xx ^ , 

9zz gzz 

2 2 

C(z) = -2^ 2 e-^ 2 %r{z)d z r{z) , M(z) = ^J^—e'^^r 2 - l)g 2 x J z . (4.32) 
gJz 27Ta x 

Notice that B(z)/A(z) takes the same value as for the vector and axial excitations, which means that 

the definition of the w-radial coordinate is the same as in those cases. The expression built from B(z), 

C(z) and M(z) which enters the Schrodinger potential takes a remarkably simple form: 

W) { M(z) ~ \ d - c(z) ) = 4 < 433 > 

We will find the UV and IR limiting behaviour of the associated Schrodinger potential. At small 
z fa u, we find H ~ /u? and one can immediately compute from ( p.7| ) the UV leading term to be 
V(u) = -£-5. Similarly, the term that dominates for large u is quadratic f^-. Thus, we have found 
that the UV and IR asymptotics are the same as for the vector case and, as a first approximation, 
we can use the soft wall equation ( [4.151) . Thus, we can again match the asymptotic behaviour of 
the current- current correlator to the perturbative result. For that, we make an argument similar to 
[ |32| assume that the qq operator is dual to the tachyon rescaled by some constant 0, such that the 
boundary coupling is, schematically, J (3{t/ z)qq. This is what we anticipated in the relation between 
Ci and the quark mass ( |3.26| ). Then, matching the large q 2 result Tls{q 2 ) = ^&<? 2 logg 2 and reasoning 
as in the the vector case, we find: 

(2ira')]CR 3 \ N c 



P 2 8tt 2 



(4.34) 



In figure [12), we depict the associated Schrodinger potential for different values of C\. Comparing 
figure [12] to figure ||, one can check that the potentials for the scalars are above those of the vectors. 
Thus, for equal excitation number, scalar mesons are typically heavier than vectors in the present 
model. 

The relation of the lowest scalar meson masses to C\ follows 

z A m { p w 2.47 + 0.683c! , z A mf w 3.73 + 0.488ci , z A mf w 4.41 + 0.507ci , (ci < 1) 
z A mf w 4.99 + 0.519ci , z A mf w 5.50 + 0.536c! , z A mf « 5.98 + 0.543ci . (4.35) 

We point out that the scalar meson masses do not depend on the parameter k. 
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Figure 12: The Schrodinger potentials associated to the scalar excitation for c\ = 0,0.5,1,2,3,4. 
4.4 Pseudoscalar mesons 

We now focus on the pseudoscalar mesons. With respect to the previous modes, there is an extra 
complication because the physical modes are a combination of two bulk fields 6 and A\. However, 
we will see that it is possible to find a combination of the fields for which one obtains a standard 
Sturm-Liouville problem. 

The quadratic action reads: 



S = -(2W) 2 /C / d^xdze-^ 2 



4 9xx9zz 

9v 



2na' 



(4.36) 



Defining (p(z) and as in (|4.2|) and Fourier transforming V(x IJ '), we can write the equations of 



motion for Ajx and 9 as: 



1 22 
-d z (e-^ 2 g xx fJdMz)) ~ k^-(<p(z) -&{z)) = 



e ^ T gl z 



/i 2 r 2 



(4.37) 



(4.38) 



These two equations can be combined into one by solving (|4.37|) for •§ and inserting this into ( |4.38| ). 

1 



e^ T T*gljJd, 



e ^ 2T2 r 2 g xx gzz 



-d z V{z) 



ll 2 T 2 



- k^^(z) - q 2 ^(z) = 



(4.39) 



where we have defined: 

ij p (z) = -e-^ 2r2 g xx fJdMz) (4-40) 

and we have used the definition of k in (|4.26|) . We can transform equation ( |4.39| ) to a Schrodinger 
form following appendix [D]. Comparing ( f4.39|) to ( p.2|) (and inserting C(z) = 0) , we find: 



A(z) 



l^r 2 -2 -l~-2 

e T 9 xx 9zz , 



B(z) = e^ 2 r- 2 g- 2 ~gl, 



M{z) 



//V 



(4.41) 
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Notice that the value of B/A coincides with those for the rest of modes and, therefore, the Schrodinger 
coordinate u is the same for all the possible excitations. Let us compute the IR (large u) leading 
behaviour of such a Schrodinger potential. It turns out to be Vjr(u) ~ ^(H- 4^)w 2 , as for the axials. 
The coefficient of this quadratic term is what controls the slope of the Regge trajectories for highly 
excited mesons. Thus, the outcome of the present model in this respect is that vectors and scalars 
have the same Regge slope, whereas the slopes for axials and pseudoscalars coincide and are larger 
than the vector one. 

An important observation is that the natural normalization condition is not flD.4|) but, looking 
for the kinetic term of the pseudoscalar field in ( |4.36|) , we obtain: 



(2ttc/) 2 K / dze 



V- 2 



A QxxQzZ 

9v 



(d z ip n (z)) 2 + 



4A 
2rca 



2 ~ 

QxxQzz 



(A 



1 

2 



(4.42) 



Rewriting this expression in terms of ip p , we find: 



1 

2 ~ 
(2na') 2 )C 
4 



(2W) 2 /C 



9 4 v 



du 



ZA 



dze^ 2r2 ( 



~gU p n {zf + 



2W 
4A^p 



g xx T 2 a(u) 2 + 



2W 



e^^gxxT 2 



d u [e"^ T gl x ra{u, 



(4.43) 



where in the last line we have changed to the Schrodinger variables following the conventions of 
appendix [B|. There are some subtleties related to the UV behaviour of ( |4.43| ) which are worth 



13 



explaining. Since the leading UV behaviour of our model is the same as in the hard wall [12 
or soft wall [|14| models, the following arguments are analogous in all these cases. However, we are 
unaware of any reference where the discussion below is explicitly shown. It turns out that this UV 



behaviour is qualitatively different for massless [c\ 
will study both cases separately below. 



0) or massive (ci ~ m q > 0) quarks. We 



4.4.1 The m q = case 

We will now study the qualitative properties of the physical spectrum for m q = 0. In this case, the 



t ~ u 3 near the UV and therefore H = (AB) * 



u 



which implies that 
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vw M = + . . . 



(4.44) 



The first correction represented by the dots comes at order u 3 . One can then find the UV expansion 
for a(u) that solves ( p.6| ) in terms of the two integration constants which we denote k\, k-i as a(u) = 
k\u~i + ^kxvn^vh — ^kxm^vh logii + A^i^ + . . . We now want to insert this in the last line of (|4.43[) 
and check whether the integral converges near u — 0. The first term is always convergent so we focus 
on the second term which behaves as J Q duu~ 3 [d u (u2a(u))] 2 . Therefore, for m n = 0, this mode is 
UV-normalizable irrespective of the values of k\ and k 2 . Thus, one can always tune kijk\ in order 
to find a solution that is well-behaved in the IR. This means that for m q = there is always a 
normalizable solution with m n = 0, which corresponds to the expected massless Goldstone boson. 
On the other hand, if m n ^ 0, one has to impose k\ = in order to have UV-normalizability. Then, 
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as in a standard Sturm-Liouville problem, there will be a discrete set of massive modes, where m n is 
selected by matching the normalizable UV and IR behaviours. 

In summary, for m q = the UV structure of the Schrodinger potential and normalizability 
condition ensures the existence of a massless Goldstone boson together with a discrete tower of 
massive excitations, as expected. 



4.4.2 The m q ^ case 

Near the UV, we now have r ~ u and therefore 



(AB) 4 ~ «2, and 



Vuviu) 



Au 2 



+ 



(4.45) 



where the first correction in the dots is O(u ). We can find again the UV solution in terms of two 

1 15 

integration constants a{u) = kyu^ log u + k 2 u^ + 0{u^). Now, requiring that the last term of (|4.43|) is 
UV-finite requires setting k\ = 0. Again, one has a Sturm-Liouville problem with a discrete spectrum. 
Figure |3| depicts a few Schrodinger potentials for the pseudoscalar mode, for different values of 




Figure 13: On the left, we plot the associated Schrodinger potential for small quark mass, concretely 
ci = 0,0.001,0.005,0.01. Even if the UV behaviour is completely different for c\ = 0, we see that for small 
ci 3 the potential looks very similar to the massless case except precisely around u = 0. On the right, we 
plot the same but with larger values of c±, in particular c\ = 0, 0.5, 1, 2, 3. All plots have been done taking 



k 
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As for the axial case, increasing k amounts to pushing up the IR part of the potential. 



Let us now look at the lowest-lying excitation when m q is small. This should be a pseudo- 
Goldstone boson with its mass given by the Gell-Mann-Oakes-Renner relation. One indeed can find 
this following the argument of [12]: for zero quark mass, the q 2 = solution of ( [4.37 ), ( 4.38 ) is given by 



d, tp(z) 



1, where we have defined ip J 



as the solution at zero momentum of 



([4.25|) with boundary condition Vy=o = 1- Consequently, regarding (|4.30|) , z l d z ip q \ =Q \ z =o = 
We now find perturbatively the small m 2 solution just by integrating fl4.38j ) and we obtain: 



6tH f2 



jj?kr : 



; d z ip a 2 =0 dz 



(4.46) 
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Using that this integral is dominated by the small z region and taking into account J ^ ~ J , g £ ( 



C32: 3 ) 2 

2^-, we can substitute the relations between ci,c 3 and m q , (qq) (|3.26|) , ( |3.27|) together with ( p^8| 
flOq) and ( g2|) to find the GOR relation [|tj 



-Am q (qq) = m\fl 



(4.47) 



We have obtained this expression by making a series of approximations. However, we can crosscheck it 
with the values for the mass obtained by the standard numerical computation, see figure |T^. Q4.47| ) is 
very accurate for small masses. When going to larger masses (up to c\ ~ 1), we can fit the mass of the 
lowest lying pseudoscalar to ^Ja m q + bm 2 q . We include here the masses of the first six pseudoscalar 
modes in terms of c\ 



(1) 

z\m p ~ 



3.53c? + 6.33ci 



(2) 

,2a my 



,( 4 ) 



2.91 + 1.40ci , 
z A my « 5.04 + 1.21c x , z A m K p> « 5.87 + 1. \lc x , 

where we have also set k = ^4 for this calculation. 



(3) 

z\my 



(6) 

z\my 



w 4.07+ 1.27ci , (ci < 1) 
6.62 + 1.15ci. (4.48) 
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Figure 14: We plot the mass of the lowest lying pseudoscalar as a function of c\ (namely, the quark 
mass). On the left, we crosscheck the GOR relation (where for {qq), f n we have introduced the result found 
numerically in the chiral limit) to some points computed numerically. The approximation is very good for 
small ci, and deviations start seeing visible around c\ = 0.03. On the right, we have fitted b in the expression 
+ bm q and checked that the fit is rather good up to c\ = 1. the parameter k has been 



-im q {qq) 
F 
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taken to be ^| in these plots. 



4.5 Mesonic excitations: a brief phenomenological analysis 

We now make a phenomenological analysis of our model by comparing our results for the spectrum 
and the decay constants of light unflavored mesons to their experimental values. An extensive study 
of meson spectrum appeared in [15|, without including the decay constants. We will fit the three 
parameters of the model, z\, c\ and k, using mesons with isospin 1 and J PC = 1 , 1 ++ , 0~ + , ++ . 
Since z A ~ ^qcd an d c i ~ m q-, h turns out that there is a single phenomenological parameter k, 
apart from those inherent of QCD physics. 
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The experimental values of the meson masses which are used are quoted by [41 . 
We fit the three parameters of our model to the masses of the light mesons which appear in table 
]T) and the decay constants appearing in ||. To make the fit we minimize the rms error 



where n is the number of the observables minus the number of the fitted parameters, n 
values of the parameters minimizing e rms read 

1 8 

z7 l = 549 MeV , c u z A = 0.0094 , k = — (4.50) 

7T 

The rms error then is e rms = 14.5% and the comparison between the experimental and model values 
appears m table d), for the masses and in table @, for the decay constants. 



jCP 


Meson 


Measured (MeV) 


Model (MeV) 


100\6O\/O 


1— 


p(770) 


775 


800 


3.2% 


p(1450) 


1465 


1449 


1.1% 


1++ 


ai(1260) 


1230 


1135 


7.8% 


o-+ 




135.0 


134.2 


0.5% 


tt(1300) 


1300 


1603 


23.2% 


0++ 


a (1450) 


1474 


1360 


7.7% 



(4.49) 
= 9-3. The 



Table 1: The results of the model and the experimental values for light unflavored meson masses. 



jCP 


Meson 


Measured (MeV) 


Model (MeV) 


ioo|<yo|/o 


1— 


p(770) 


216 


190 


12% 


1++ 


ai(1260) 


216 


228.5 


5.8% 


o-+ 




127 


101.3 


20.2% 



Table 2: A comparison of the results to the experimental values for the decay constants of light unflavored 
mesons. 



5. Meson melting in the deconfined phase 



We briefly discuss in this section the fate of the mesonic modes when the gauge theory undergoes a 
deconfining phase transition | 42]| , namely when we use the background of equation ( |2.8|) . The first 
observation is that, as we saw in section |3.2j , the tachyon cannot diverge at any point in this case 
and, therefore, the brane reaches the horizon, and we only have "black hole embeddings" , in analogy 
with the terminology introduced in |43|. This means that there is no discrete spectrum above the 



deconfining phase transition. When we are considering small quark masses, this is perfectly realistic. 
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However, in the real world, charmonium and bottomonium do survive the QCD phase transition. 
We want to study this problem in the present model, and therefore we will compute the spectral 
functions at different values of m q /T. In particular, we will focus in the vector excitation. 

We start by discussing the associated Schrodinger potential for the vector excitation in the de- 
confined background, at zero momentum. The expressions in (|4.7| ) are modified to: 

A{z) = e~-^ 2 gLAfJ , B{z) = e'^gLg^-gl , (5-1) 

where one should remember that now refers to the metric (|2.8|) . Notice that B/A diverges at 
z = zt as a single pole, such that J ^ B j Adz diverges and the horizon z = z T corresponds to u = oo 
in the Schrodinger coordinate. V(u) is exponentially decreasing for large u. For completeness, we 
write in appendix ^| the functions determining the potential for the rest of modes. In figure [15], we 
show several examples of potentials computed numerically, for different values of C\ ~ m q ZT ~ m q /T. 
In the second and third plots, we also compare it to the potentials in the confined phase for the same 
value of c\. Namely, we show how the potentials for the vector excitations are modified at the phase 
transition. They share the same UV behaviour, but are drastically modified in the IR due to the 
different behaviour of the tachyon and the metric. 




Figure 15: The Schrodinger potentials associated to the vector excitation in the deconfined phase, at zero 
momentum, for different values of c\ ~ m q /T. In the first plot c\ = 0.01, 1, 2, 3, 4. The second and third plot 
(respectively c\ = 1,2) make a comparison with the potentials in the confined phase for the same values of 



In [46]], it was shown that a step-like potential gives quasi-particle behaviour. Moreover, if the 
position of the step coincides with that of a barrier in the confined phase, the quasi-particle mass is 
related to the mass of the meson before the phase transition. The third plot seems to point along 
that line. However, in the present model the potentials for the deconfined phase are never step-like 
enough, nor present bumps, and thus do not create sharp peaks in the spectral function, as we will 
see below. 

Once we have the potentials, it is straightforward to compute the spectral function from the 
retarded correlator which is computed following the prescription of [47||. In practice, what one has 
to do is to impose ingoing boundary conditions at the horizon (a(u) ~ e lulu near u = oo) 8 and find 



8 Even if we use the notation with u since it is better for illustrative purposes, we have found easier to perform the 
numerical computations in the z-coordinate. 
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the behaviour of the wavefunction near the boundary, namely compute &i, 62 matching the numerical 
result to the UV-expansion ( f4.12|) . Then, the spectral function is given by: 



-Im Gr(uj) 



(5.2) 



where we have replaced the parameters of our model by J^, using ([4.20 ). We show in figure 16 



examples of this computation. We have plotted j^^t- in terms of u for various values of c\. 



some 




Figure 16: The spectral function (divided by u' 2 ) in units of M± for c\ = 0.01, 1, 2, 3, 4, from top to bottom 



In ||44| , a bottom-up holographic model was built in order to discuss the physics of the J/^f above 
the phase transition (see also for related work). In the model of |P| , there is also no discrete 
spectrum but, for low enough temperature, by engineering the potential, the J/\l/ shows up as a 
peak in the spectral function (or, alternatively, as a small negative imaginary part of the associated 
quasinormal frequency). Namely, it has quasi-particle behaviour. 

We now discuss the results comparing to 



The authors of [44] showed how the Schrodinger 



potential should look like to find some qualitative physical properties of the J/\l/ meson. In our 
model, the potentials are found dynamically, and look similar to except for the important fact 
that they fail to have the dip of ||44|| . As a consequence, the J/\l/ decay constant is too low (see figure 



t) and a clear peak in the spectral function is absent. One may compare our figure [16] to figure 3 in 



It does not seem possible to introduce the dip feature in our model without inserting at least a 
new scale in the problem. Notice that the Aqcd does not play any role in the deconfined case for the 
present model, apart from setting the value of T at which the transition takes place. In other words, 
the z\ does not appear at all in the deconfined background. We speculate that one should modify the 
background by including some dependence on Aqcd in the deconfined phase in order to improve the 
model in these respects. It would be interesting to investigate whether analyzing the tachyon action 
in the less simple but well motivated backgrounds of "improved holographic QCD" models [[JIJ might 
ameliorate the issue along this path. 
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6. Discussion 



6.1 Summary and general comments on the model 



In this paper, we have analyzed in detail several issues of the holographic model presented in |15 



The main ingredient is to describe the open string (meson) physics by using tachyon condensation, 



captured by Sen's action |fL7l , as first advocated in ||16|| . The "tachyon", namely the lowest lying bifun- 
damental scalar of the brane-antibrane system, is dual to the quark mass scalar bilinear operator and 
its condensation corresponds to chiral symmetry breaking. In order to make explicit computations, 
it is necessary to give an explicit form of the tachyon potential and to choose a curved gravitational 
background. We have considered extremely simple possibilities for both, see ( |3.4| ) and ( |2.2| ). 

As we have already remarked, the model is inspired by string theory but is phenomenological. It 
does not provide a well controlled approximation to string theory in any limit. In this sense, it should 
be considered a bottom-up model. Nevertheless, the main point we want to make in this work is that 
it can be very useful to incorporate top-down derived ingredients as (|3.1|) into bottom-up models. We 



will make a comparison between qualitative properties of our model and the soft wall model and 
modifications thereof. Successes of |14| include a reasonable qualitative and quantitative matching of 
properties of the measured mesonic states, including the Gell-Mann-Oakes-Renner relation and the 
Regge trajectories. These are also found in our model, together with the following extra appealing 
properties, which are automatic in our general set-up: 

• The model incorporates confinement in the sense that the quark-antiquark potential computed 
with the usual AdS/CFT prescription [[|9j confines. Moreover, magnetic quarks are screened. 



The background solution stems from a gravitational action, that allows, for instance, to compute 
thermodynamical quantities. All of this are properties associated to the background geometry 
and were already discussed in [[HI] , 

The string theory nature of the bulk fields dual to the quark bilinear currents is readily identified: 
they are low-lying modes living in a brane-antibrane pair. 

Chiral symmetry breaking is realized dynamically and consistently, because of the tachyon 



dynamics. See [48] for discussion and possible solutions in the soft- wall model context. 



In the present model, the mass of the p-meson grows with increasing quark mass, or, more 
physically, with increasing pion mass. This welcome physical feature is absent in the soft wall 
model, fl4| . It occurs here because the tachyon potential multiplies the full action and in 



particular the kinetic terms for the gauge fields, which therefore couple to the chiral symmetry 



breaking vev. In our previous work [15], we exploited this fact in order to fit the strange-strange 



mesons together with the light-light mesons, with rather successful results. In |50| , the authors 



added the strange quark mass to the hard wall model and computed the dependence of vector 
masses on the quark mass. In that case however, this dependence of the vector masses originated 
only from the non-abelian structure and therefore misses at least part of the physics 9 . 



9 On the other hand, the quark mass dependence of the p-meson can be seen in different top-down models, see |HJ 
for a recent work in the context of the Klebanov-Strassler model. 
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• The soft wall requires assuming a quadratic dilaton in the closed string theory background. It 
has been shown that such a quadratic dilaton behaviour can never be derived from a gravi- 
tational action while keeping the geometry to be that of AdS. 10 . That the background is not 
found as a solution is a shortcoming if for instance one wants to study the thermodynamics of 
the underlying glue theory. The thermodynamics of the soft wall model is therefore ill-defined. 
In the present model, we found the background as a solution of a two-derivative approximation 
to non-critical string theory, see section Q In order to obtain Regge behaviour, we also needed 
a further assumption: that the tachyon potential is asymptotically gaussian. However, this is 
rather natural since this potential has appeared in the literature, for instance |26|, . Still, we 
warn the reader that the formalism of |26|, cannot be directly and controllably connected to 
the present setup. 



Considering that the dynamics is controlled by a tachyon world-volume action automatically 
provides the model with a WZ term of the form given in |52, 53, 27\. We have not discussed 
this term at all in this work, but in Jl6| it was shown that properties like discrete symmetries 
(parity and charge conjugation) and anomalies are, in general, correctly described by analyzing 
this term. 



In summary, we regard our model as being in the general framework of [jnfl , but with several 



qualitative improvements due to the dynamics built in the action (|3.1|). Moreover, our starting 



assumptions are rather simple and well motivated from a top-down perspective, such that the ad 
hoc input is scarce. It is also encouraging to find that quantitative fits to experimental data are 



reasonably good, see [J5] and section [4.5| , but those are not the main aim of the present work. 

We have also found several aspects in which the model does not capture features which are known 
from perturbative QCD. 

The main issue is the leading large q 2 contribution to IIy(g 2 ) — n^g 2 ) as will be discussed in 
section |0|. We have also seen that our model does not work that well for large quark masses 11 as 
it would grossly underestimate the decay constants for charmonium and no clear peaks are observed 
in the corresponding spectral functions in the deconfined phase (see the discussion of section fty. It 
would be interesting to know whether mild corrections of the model could ameliorate these issues or 
whether these are unsurmountable differences of models of this class to actual QCD physics. 

The behaviour of for small m q (figure |2|) has been investigated in earlier classic works [54 . 
However, the leading IR divergence-free correction, is of order 1/N 2 at large N, originating from a 
pion loop, |j55l. The leading-N corrections come from the four-derivative terms and are dominated by 



the scalar meson contribution. The behavior is qualitatively similar to what we find. For large m q we 
have not been able to calculate the asymptotic behavior but generically speaking we do not expect 



10 This was shown in the second reference of [||. In Q| such behavior can be implemented for glue, but the metric 
changes appropriately, an important ingredient for implementing confinement in the glue sector. 

11 Perhaps this is not surprising since for heavy quarkonium perturbative methods and in particular non-relativistic 



QCD (see |56 for reviews) are accurate and it may be naive to expect that a dual theory can be a good approximation 



to the physics. 
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it to necessarily match that of QCD, as discussed in [54]. The reason is that the UV asymptotics of 
the bulk gravity solution are not necessarily the same as in QCD. 



6.2 Comments on effective actions for the open string tachyon 



The notion that a scalar bifundamental in a brane-antibrane system should be the holographic dual 
of QCD-like chiral symmetry breaking is rather simple and robust, see for instance |TF| ^8|. |60| . 
What is not obvious, however, is which effective action is best in order to describe a brane- 



antibrane system if curved spacetime. We have used the simple proposal of [[17], but one should keep 
in mind that other alternatives might also be useful. We provide here a short guide to the literature 
on the issue. 

Garousi and collaborators, starting from the early work [pl}| , have tried to use explicit string 
theory computations in order to constrain the tachyon generalization of the DBI and WZ actions 
[p^| , see also f6"3fl . In |p8| , an action for a Dp- Dp system based on a particular symmetriced trace 
prescription was proposed 12 . There are subtle differences between the proposal of p8f and the one by 



Sen [17|, which may have dramatic consequences. In fact, we have checked that using the symmetrized 



trace action of [[28| for our model, one still finds Regge trajectories for vector and axial mesons but 
the slope for the axials changes, see appendix ||. The study of other physical properties inferred from 



the action of [61] is beyond the scope of the present work. 



In ||64|| , it was discussed how to take into account the brane-antibrane distance in the string action. 
This is important for holographic duals, specially if one wants to insist in top-down approaches. A 
reason is that, in the weakly coupled picture, if one has brane and antibrane on top of each other in flat 
space, the tachyon would create a real instability (it cannot be compensated by the AdS curvature). 
Therefore one should think of separated branes as in the Sakai-Sugimoto model [O. With this in 



mind, generalizations of [11], based on the action proposed in [54] where constructed in | 65fl . 

In a beautiful recent paper |66| , building on the work |67fl , Niarchos proposed a different way of 
building the effective action which should better capture the physics of a separated brane-antibrane 
system. This may be useful in holographic modeling and in particular to improve the Sakai-Sugimoto 
model. We also refer the reader to for a more exhaustive overview of the literature of tachyon 
effective actions. 



6.3 On the OPEs and the slope of the Regge trajectories 

There has been some debate in the literature on whether the large Euclidean momentum behaviour of 
two-point correlators can be used to constrain the behaviour of the QCD spectrum of excited mesons. 
The main point is to compare infinite sums like (|4.21|) (or, more precisely, differences of such sums: 
vector minus axial or scalar minus pseudoscalar) to the large q 2 behaviour expected from the operator 



12 In [pSfl , a trace is needed even for a single brane-antibrane pair since the degrees of freedom are 2x2 matrices. In 
order to non-abelianize the flavor group in our case, one should also deal with the problem of how to implement traces 



on a non-abelian generalization of Sen's action (3^). Investigating the physical consequences on the dual theory of 
this non-abelianization and of different proposals for the effective action p2| would be very interesting, but is beyond 
the scope of the present work. 
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product expansion (OPE). In particular, there is the question of whether different Regge slopes in 
the vector and axial (or scalar and pseudoscalar) channels 



(m v J A ? ~ 



Ay^n + const for large n 



(6.1) 



are consistent with the OPEs. Notice that this is a theoretical question, irrespective of the exper- 
imental observation of the spectra. Let us give a brief and incomplete overview about the debate 
regarding this issue. For instance, in [jSfl, a model with Ay 7^ was put forward. Later, in |6{J, it 



was claimed that this model was inconsistent with the OPEs, but the arguments of |)9] were called 
into question in | 70fl , due to subtleties in the regularization of infinite sums like ( |4.21| ). More recently, 



works like J7I], |72|, [7J| claim that the Regge slopes should be equal whereas the opposite conclusion 



was reached in [[73], [75], |76 |. 

We have found above that in our model, there are different asymptotic Regge slopes A^ > Ay. 
However, the coefficients of the leading logarithms in the large q 2 correlator for vectors and axials 
coincide, consistently. In order to illustrate this fact, let us remember that, asymptotically in the 
UV, our model resembles the soft wall of |14[], see equation ( |4.15| ). In the soft wall model, the Regge 
slope is controlled by the constant c of ( |4.15| ), but the quotient F 2 /A 2 is independent of c f77 |. 
This quotient is indeed what controls the coefficient of the leading logarithm |FJJ. We thus have 
lim^^.^ (LV(g 2 ) — YlA^q 2 )) = together with different Regge slopes. However, this is not enough 
to comply with the OPEs. The leading contribution to (Uv{q 2 ) — Ha{q 2 )) at large q 2 should be of 
order q~ 4 because QCD does not have dimensionful quantities that allow to rewrite for instance a q~ 2 
term. We can resort to numerics to compute (Hy(g 2 ) — n^g 2 )) in our model and the result does not 
comply with the q~ A expectation. This fails in the axial channel as shown in appendix |F|. 

The obvious guess is that, since our holographic model is clearly not exactly QCD, it includes 
operators or condensates which are absent in QCD, modifying the subleading pieces of the correlators. 



The same kind of problem is present for any holographic model we are aware of, see ^Tt§ for a 
discussion concerning the soft wall 13 . However, this seems to us more of a technical problem that may 
be resolved by finding the appropriate potentials than a general obstruction to this class of models. 
Settling these issues requires further work. 

6.4 Outlook 

As we have discussed, our simple model is quite successful in describing many features of QCD. A 
lot of effort has been devoted in bottom-up models to estimate other QCD related properties as for 
example form factors, see for instance [ 78| . To reproduce such computations in the present setting 
and compare the results is an interesting problem that we leave open for the future. We have not 
studied non-trivial baryon number or chemical potentials, which would clearly be worthy extensions 
of the model. 

There are some aspects of the present that would be worth improving, like the physics associated 
to heavy quarks (compared to the QCD scale). We have just explored the result of working with 
Sen's action in the simplest confining holographic background available in the literature 



20, 21 



13 



A. P. thanks O. Cata for a discussion on this subject. 
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Therefore, it is still left to understand the consequences of implementing the tachyon action in different 
backgrounds, as for instance those which go under the name of improved holographic QCD |35| or 
modifications thereof. 

It could also be interesting to try to introduce the quarks beyond the quenched approximation 
and therefore compute the backreaction of the tachyon action for the fundamental fields on the 
gravity background. For a review of unquenched flavor in critical (ten-dimensional) string theory 
backgrounds, see [[7jJ. 

Of course, it would be worth to provide a non-abelian generalization of this model, for which one 
should provide a technical prescription on how to take traces in the action. Another line of obvious 
interest would be to use the model for baryons. Since quark masses play a more dynamical role than 
in other bottom-up approaches, this could be interesting for the physics of the sigma-term. 

Finally, and most importantly, we would like to point out that using an effective open string action 
like ( |3.1|) in the framework of holography can well have interesting applications beyond the realm of 
strong interactions. For instance, many bottom-up technicolor models have appeared in the literature, 
see [BUI for a review. It is a very interesting question to understand whether chiral symmetry breaking 
controlled by an action like fl3.1p may offer new dynamical possibilities for the modeling of electroweak 
symmetry breaking. On the other hand, in the last years, many phenomenological models have been 
constructed in order to address some issues of superconductors and other condensed matter systems, 



for a review see |pl |. Again, we would like to remark that (|3.1| ) could hopefully lead to interesting 



new dynamics in different set-ups. 
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APPENDIX 



A. Book-keeping summary of the parameters of the model 

We summarize here the parameters of the our model. There are two parameters coming from the 
background, the AdS radius R and the position of the cigar tip za- 

The action of the flavor brane-antibrane pair also includes a' and two more parameters gy and 
A which are related to the normalization of the vector and tachyon fields respectively. The tachyon 
potential also includes two constants, /C which is an overall factor in front of the action and /i 2 . It 
should also noticed that \i can be absorbed in r(z) by redefining f — > fir. Then this parameter 
disappears from all equations. We used /i 2 = it, for the numerics through our analysis, but this does 
not affect any physical results of our model. 

Another parameter which exists in the model is c\ which appears in the UV asymptotic of the 



tachyon expectation value ( |3.7|) . C\ is proportional to the quark mass, with a proportionality constant 



f3 which was introduced in (|3.26|) , however (5 does not appear in the equations for the spectrum or 
the decay constants, so its value is not relevant for the model predictions. 

In total we have the following parameters R, za, a', gv, A, fC and c\. R, a' and A are related by 
equation (|3.8| ), which relates the tachyon mass to the dimension of its dual operator. Then, we relate 
gv and A to the number of colors N c in QCD by matching the results of the vector and scalar two 
point functions as calculated in bulk on the one hand and in QCD on the other hand. The results 
are given in ( [4.20| ), ( [4.34j ) , which relate gy to A. Hence, finally the spectrum and the decay constants 



depend on z^, C\ and a combination of R, gy and a' which was named k and is given by, Q4.26D 



B. Analysis of singularities in the tachyon differential equation 
B.l Confining Background 

In this appendix, we will investigate the existence of singular solutions of the tachyon equation of 
motion (|3.9| ). It was argued already that r can diverge only at the tip of the cigar. Therefore what is 
left to investigate is solutions where r'{z) diverges at a point z G [0, z\], but r{z) remains finite at the 
same z . We call these solutions "spurious". Taking into account t'(z) ^> t(z) in the neighborhood of 



zq, the leading terms of Eq^^) are the first and the second one. Hence, the leading order equation 



is (we set /i = 1 in the sequel as it can be absorbed in the normalization of r and z\ = 1) 

t"{z) - \zf{z)r{z) = (B.l) 
with solution in the vicinity of the divergence 

t'{z) = -±== 1 , T {z)= f -p= + r Q (B.2) 

V9(z) ^C-|^ 2 (l-fz 5 ) Jo V9(z) 
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where g(z) = C — |z 2 (l — f-2 5 ). The function g(z) has either one or three real roots. In particular, 
there are the following three cases 



1. C < : There are no roots of g(z) in the interval [0, 1], since the only one root is at Zq > 1. It 
should also pointed out that for z £ [0, 1], g(z) is negative so the solution does not exist. 

2. < C < |y : There is a single real zero at z £ [0, 1]. While, the other two real zeros lie outside 
that interval. 

3. When C = |^ the divergence happens exactly at the tip of the cigar. 

4. C > |y : Again there is no real zero in [0, Zq]. 

If g{z) has a real root zq £ [0, 1], then it follows from ( |B.2|) that t'(z) diverges at z = zo, but 
t(z) is regular there. Only in case that z is a double root of g(z), both t(z) and t'(z) diverge at the 
same point. We are particularly interested in the above case where the acceptable solution diverges 
at some point z in order to obtain the effect of chiral symmetry breaking in the dual quantum field 
theory. This is only managed if we tune the initial conditions (C here). There are two possibilities 
which lead to a double root of g{z), in the context of the above approximation 



1. C = : In that case the double root is at z = 0. Then, Eq.( p.l| ) has not real solutions, so it 
is not considered here. 



2. C = |y : The double root now is at z = 1. 

The only rigorous result of the aforementioned analysis is that "spurious" singularities are generic 
if < C < § . 

If the tachyon diverges as a power law r ~ (z — zo)~ a with a > 0, ^ ~ an d the approximation 
described above is still valid, provided zo ^ 1. But then, this is not a valid solution since the solution 
in the above approximation is of the form ( [B.2| ). This excludes such a divergence if Zq ^ 1. 

The only other option of divergence of t(z), and/or t'(z) at a point Zq £ [0, 1] is the case where 

term is dominating. The relevant equation then is 

t"(z) + r\zfr = (B.3) 

which leads to 

r \z) = C e-^ 2 (B.4) 
For t(z) ^> 1, an approximate solution to Eq.( B.3| ) is 

_L e \<z? ~Cz + ... (B.5) 

t(z) 

Therefore, t(z) diverges only if z — > oo which is not allowed in the present geometry as z £ [0, 1]. 
Hence this case is excluded. 
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From the above mentioned, we conclude that the only place where r(z) diverges is at z = 1. In 
order to find the solution of diverging t(z) we must tune the initial conditions in the UV, see section 
|3|. We also showed that "spurious" singularities in the interior of the interval [0, 1] are generic for a 
range of initial conditions. 

The existence of "spurious" singularities has been verified numerically, and it fits the asymptotics 
( p.2| ). An example of this behavior is shown in Fig. [17|. We solve numerically eq. (|3.9| ) for arbitrary 
initial conditions, meaning that the mass and the vev are not tuned according to the plot in Fig. 
|2|. In particular, we have chosen c\ = 0.1, C3 = 0.439. In this case we notice that t'(z ) ^> t{zq) 
at zq = 0.8696 < 1. The right part of Fig.flTTD includes the plot of the derivative of the numerical 
solution (red line) and the expression for r'(z) given in Eq.( |B.2| ) (dashed blue line), near zq. On the 
left part we have plotted r(z) and the asymptotic solution ( |B.2| ) (blue dots). The parameters of the 
expressions in ( P-2|) are C = 2.71758 < l^vr and r = —0.28. For those values the asymptotic solution 
QB.2j ) fits the numerical solution of the full equation near z . 





Figure 17: An example of a solution of type which diverges at zq = 0.8696 < 1 was found numerically. The 
numerical solution and its derivative are compared to the asymptotic solution 



B.2 Deconfined Background 

We now look for singular solutions of the tachyon equation of motion in the deconfined background, 
Eq.( 3.17j) . Considering again that t'(z) ^> t(z) at the vicinity of Zq, with z$ being the point where 
t'(zo) —> 00, the leading terms of Eq. (|3.17|) near z are 

X»W4/W (-1 + ||>)a^0, (B.6, 
where f(z) = 1 — The solution reads 

t(z) = f dZ = + r (B.7) 

3 Z 5 



C-\z- (l- 



28 zg 
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The function g(z) = C — ^1 — J|f>r^ has a maximum in z = and no other extrema in the 
interval [0, zj\. If this function is zero at some point Zq, then t'(zo) is infinite and t(zq) is finite except 
if g(zo) = has a double root. Similarly to the confining background case, there are three choices 

1. C < : g(z) is not zero in the interval [0, zt\- 

2. < C < |f Zt : d( z ) nas one rea ^ ro °t m the interval [0, zt] which is not a double root. When, 
C — |f the root of g(z) is at z = 1. 

3. C > ff : There is no root of g(z) in [0, zj\. 

So, for a suitable range of the initial conditions (C) we may have a solution with diverging t'(z) and 
finite t(z) at some point z G [0, 1]. A double root is possible to be found in case of C = 0, and it is 
at z — 0. 

The discussion about the case where the term r' 2 r of Eq. fl3.17|) is leading, remains the same as the 



one that follows Eq. ( |B.3| ) in the previous appendix. Therefore, in case of the deconfined background 
tachyon cannot diverge in the interval [0, 1] but "spurious" singularities of the form r'{z) 3> t(z) 
exist in general. 



C. Scheme dependence of the condensate and the constant a 

We have introduced an arbitrary constant a in (|3.24| ) associated to a counterterm that gives a finite 
contribution and therefore cannot be fixed by demanding the cancellation of divergences. This is a 
common situation in holographic renormalization and is related to scheme dependence of renormal- 
ization in the field theory side, as we now discuss. 

The gauge invariant composite operator q\qR must be defined by a subtraction in QFT. To asses 
what enters in such subtractions we must study the OPE q\ j {x)qR(y) as x — > y. The operator can then 
be defined by point splitting, subtracting divergent contributions, and then taking the limit x = y. 

Apart from the identity operator, all other operators that can appear in the OPE q\ j (x)qR(y) do 
not provide divergences. These include the operator itself : q\ j {x)qR(x) : as well higher dimension 
operators. Therefore this composite operator can be defined by normal ordering. To make this precise 
we use a Dirac basis, so that we have the real part q{x)q(y) and the imaginary part q{x)rfq{y). The 
imaginary part is finite and no subtraction is needed. For the real part we introduce a momentum 
cutoff A as we will be working in momentum space. We will therefore define 

: qq := lim 

A— >oo 

where on the left-hand side we subtracted the most general expression of scaling dimension 3. This 
should be enough in a conformally invariant theory. If the theory is asymptotically free, as in QCD, 
more subtractions are necessary, as there is one more relevant scale entering the problem, namely 
Aqcd- However in the model we consider in this paper, the physics in the UV is conformal so these 
subtractions are enough. 



A 2 

qq — aiA 3 — a^mA 2 — a^m 2 A. — c^m 3 log — - — 05m 3 

m 1 



(c.i; 
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To establish the coefficients in (|U.lj) we must require that the (perturbative, short distance) part 
of the vev is finite when we remove the cutoff. We have 



(: qq :) = lim 

A— >oo 



A 2 

(qq) — aiA 3 — 02mA 2 — a 3 m 2 A — a^m 3 log — - — a 5 m 



We calculate (in Euclidean space) 



(qq) 



d 4 p —ifi + m 



-N Tr 

' ' (2vr) 4 g 2 + m 2 



-AN r m 



d A p 



(2vr) 4 g 2 + m 2 



-2N c tt 3 m 



A2 p 2 dp 2 
p 2 + m 2 



-2N c n 3 m 



-2N c tt 3 m 



A 2 



A 2 - m 2 lo. 



2 , A 2 + m 2 



A 2 - m 2 log — + 0(A 



-2\ 



where ^3 is the volume of th unit 3-sphere. 
To renormalize we must choose 



(C.2) 



(C.3) 



a\ = 



a 2 



-a 4 



a 3 = 



(C.4) 



while 05 can be arbitrary. It is this arbitrariness that is reflected in the coefficient a in the holographic 
renormalization setup in (|3.27|) . 14 

In some cases, this scheme dependence can be fixed on physical grounds. For instance, in su- 
persymmetric cases, one can demand that the renormalized on-shell action is zero, see |82| for an 
example. We could not find any convincing prescription to fix a in the present case. A possibility 
would be to demand that the condensate vanishes for asymptotically large quark masses, but we have 
checked numerically that it is not possible. This is not surprising, since we have seen that our model 
works much better for small masses than for large ones. 



D. Appendix: converting to Schrodinger form 

In many cases, it is useful to write down the Sturm-Liouville problem which determines the spectrum 
of any given mode as a Schrodinger-like equation. Let us start by writing a generic quadratic five- 
dimensional action for a field ^/(x^, z): 

S = ~K 9 J d 4 xdz (A(z)(d z mf + B(z)^ u d^d u m + C{z)^d z m + M(z)^ 2 ) , (D.l) 

where we have allowed an arbitrary constant multiplying the action. Let us consider \1/ = e tqx il)(z) and 
define as m 2 the discrete set of values of — q 2 which satisfy the appropriate normalizability conditions 

14 In a theory like QCD, the arbitrariness involves the addition of a finite function m 3 f ^ Aqcd ) renec ^ s the 

presence of an extra mass scale. 
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of the Sturm-Liouville problem. The discrete set of solutions satisfy the equation of motion extracted 
from flPlD : 



Biz 



-d z (A(z)d z ip n (z)) + h(z)ip n (z) = m 2 n i/j n (z) 



where we have introduced: 



Kz) = ^- (m(z) - l -d z C{z) 



We can define the orthonormality condition: 



dzB(z)ip n (z)i) m (z) = 5 r , 



We now define a new radial variable w, and a rescaled field a in terms of a function H as: 



du 



'Biz) 
A(z) 



dz , a = E tp 



(A(z)B(z)y 



The Sturm-Liouville problem now takes the Schrodinger form: 

d 2 a n (u) 



du 2 



+ V(u)a n (u) = m n a n (u) 



where the Schrodinger potential is: 



V(«) = ig + A(«) 



(D.2) 
(D.3) 

(D.4) 



(D.5) 



(D.6) 



(D.7) 



Substituting ( |D.4| ) in (p.5| ), we find that in the new variables, the normalization condition is the 
canonical one: 

dua n (u)a m (u) = 5 mn (D.8) 



In order to estimate the mass of the modes with large eigenvalues, it is sometimes useful to employ 
a WKB formula: 



d(m 2 n 
dn 



2tt 



du 



-i -i 



ui \/m 2 n - V(u) 



(D.9) 



where U\ and u 2 are the classical turning points. 



E. Two-point function and sum rules 

We describe here how, typically, the bulk solution needed to compute a two-point correlator from the 
gravity side can be written in terms of an infinite sum. Physically, two-point functions at arbitrary 
momenta are expressed as a sum over the discrete set of physical states. This is not a new result 
but, however, we believe that explicitly making the discussion as shown below can be illustrative. We 
remark that the argument of this appendix is not general in the sense that we have not included in 
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the reasoning the possibility of having counterterms or other subtleties which may to be dealt with 
in a case by case basis. 

Let us start with the equation: 



Biz 



-d z (A{z)d z i>{z)) + h{z)if>{z) = -q 2 ip(z) 



(E.l) 



where we have used the notation of appendix [D]. We want to find a solution ip q (z) such that in 
the boundary ip q (0) = 1, and which is IR-normalizable. Our goal is to write ip q (z) in terms of the 
discrete infinite set of solutions ip n {z) of the Sturm-Liouville problem (|D.2|) , normalized as ( p.4|) . Let 
us momentarily change to the notation with a(u), in which the problem is converted to: 



d 2 a{u) 
du 2 



+ V(u)a(u) = —q 2 a{u) 



(E.2) 



and the discrete spectrum is a n (u) with (p.8|) as normalization and the completeness relation: 

a n (u)a n {u) = 5{u — u) 



(E.3) 



We introduce: 



a n (u)a n (u') 
q 2 + m 2 n 



(E.4) 



such that it is a Green function, 



du 2 



— V(u) — q 2 G(u,u') = 5(u — u'). Let us assume now that the 



UV boundary is at u = and that UV-normalizability implies that a n (0) = 0, such that G(0, u') = 0. 
Regarding ( p.5|) , a q = H^, such that for generic momentum the UV condition is a q (0) = H(0). The 
solution we are looking for reads: 



a q (u) = E(u) + / G(u, u')(h(u') + q 2 E(u'))du' 
Jo 



(E.5) 



We can translate this back to the original variables. After some manipulations, we get our final result: 



?Pn(z) 



m 



n JO 



m 2 n (q 2 + m 2 n ) 



(E.6) 



Two-point correlators are built from the on-shell action associated to this solution, which can be 
typically found by computing the derivatives fo i/j q (z) at the boundary. From the last term in (|E.6|) 



one can find the decay constants of the states in the spectrum, as in section \j.l.3[ The second term 
is g-independent and in fact it can be thought of as the q 2 = contribution. For the axial excitation, 
this is related to the pion decay constant whereas for the vector excitation, this term is absent since 
h(z) = in that case. 
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F. Axial current-current correlator 



The axial current-current correlator is now derived following the same procedure as in section ( [4.1.2|) . 
We are interested in the two point function in the limit of large Euclidean momenta. We expect 
that the leading term will be the same as in the vector current correlator, but the subleading term is 
different. 

We define the correlator as 

f e*** < J,(x)J„(0) >= (q\ u - q,q u )U A (q 2 ) (F.l) 

J X 



This is calculated by differentiating twice the on-shell bulk action. Integrating by parts (|4.24|) we 
find 

S A ,re 9 = ^J^/C / (e~^ 2 g xx fjA^ z)d z A»(-q, zj) (F.2) 

9v J ( 27r ) v ' z =£ 

where A^q, z) = if) A (q, z)Ao(q). Then, we insert the asymptotic solution ( |4.28|) into the action 



Sa = ^^-KR j -^Al{q)A»{-q) (26,(9) + (q 2 + k^c\){\ + loge)) (F.3) 

where we have set b\ — 1. The last term is cancelled by the corresponding counterterm from fl4.4|) , 
so after differentiating twice with respect to Ao(q) we find the final answer 

9v T 

We now compute 62(9) for large q 2 , similarly with section (|4.1.2| ). We convert ( |4.25|) to Schrodinger 
form. Then, the new variable u reads 

u ~ z , u — —j=[it (F-5) 
V3 

in the UV and IR respectively. We calculate the asymptotic behavior of the Schrodinger potential in 
the UV and IR 

V uv ~ £j + kfi 2 cj , V IR ~c 2 u 2 (F.6) 

where c 2 = -x (§ + 3k). Adding these two contributions we finally find the equation of motion the 
axial modes 

-d 2 u a + (J^ + kii 2 c\ + cVj a + q 2 a = (F.7) 
where a(u) ~ u"^ip A {u) near z — 0. Its general solution reads 
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a („) = fefXp C £-t^M, 0, + fe^U- 1 C £±*^, 0, «" ) , (F.8) 
/it V 4c / \/it \ 4c 



we set &2 — 0, since the generalized Laguerre polynomial is going to infinity in the IR. k% is found 
such that lirn^o ip (q, z) — 1 



4c V 4c / V 4c 

By expanding the solution for large momenta we find 



6 2 1 2 1 ^ 2 c 2 (2 7 -log4 + logg 2 ) 3fc 2 // 4 cf - 8c 

q 2 ^oo q l 4 4 24g 4 

and eventually 



/Ci?(2W) 2 / 2 , i „ o ^ ,Vc 2 (2 7 -log4 + logg 2 ) , 3P/i 4 c 4 - 8c 2 



IW) = lQ Sg -(l + log4-2 7 )+4 ^ 1V ' f ^+ j 

(F.ll) 

We notice that for c\ = the result coincides with the vector current two point function. The 1/g 2 
term, which is absent in the QCD result, comes from the mass term of the axial field (~ rA^#), 
see (fO§. 



G. Excitation equations in the deconfined phase 

We assemble here the the Shrodinger functions, as defined in Appendix B, for the equations of motion 
of field excitation modes in the deconfined background fl2.8|) , at vanishing spatial momentum. These 
modes satisfy an equation of the form ( p.2|) but there is no discrete spectrum. 

In case of the vector excitations, we have already mentioned the functions giving rise to the 
Schrodinger potential approach in ( |5.1|) . Then, the variable u reads 



u=[\l^§-dz, (G.l) 
Jo V 9tt{z) 

which remains the same for all different excitations. 

For axial- vector mesons A(z) and B(z) are the same as for the vectors (|5.1|) , but now we also 
have: 

h(z) = ^ = k^f T {z). (G.2) 
The Schrodinger functions appearing in the equations of motion of scalar excitation modes are 

A(z) = e-^glj^ , B{z) = e -^gLg^%, 

g Z z 9zz 

C{z) = -2^e-^ 2g -^r{z)d z r{z) , M{z) = ^-e'^^r 2 - l)glj t gk . (G.3) 

gll 27r « A 
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The function M(z), B(z) and C(z) combine to give a quite simple expression for the h(z) defined in 

my- 

h(z) = ~f T {z) (G.4) 



z 2 ' 



Finally, for the pseudoscalars we have 



A(z) = e^ 2 r- 2 gjgk^ , B{z) = e^V 2 ^!* , h(z) = = k^- f T (z) . (G.5) 



M{z) _ uV 
B(z) ~ z 2 



H. The action with the symmetric trace and Regge slopes 



In [28], Garousi proposed an effective action for the brane-antibrane system which has subtle difference 



with respect to Sen's one [|17|, which we have used in this work. One may wonder what would be 



the physical consequences of using such an action in our model. We focus in this appendix on the 
behaviour of the spectra of highly excited vectors and axial vectors. The equations for the vectors 
are not modified with respect to the main text. The equations for the axials are different. It turns 
out that they still obey a Regge law m 2 cx n for large excitation number n but with different slope 
compared to the main text. This slope is still larger than the one for vectors. We will not deal in the 
present work with other possible phenomenological implication of this different tachyon action. 
Garousi's action reads 15 : 

S = -STr J d 4 xdz e~ f2 \]- det(g MN + F MN + D M TD N f) (H.l) 

where hatted symbols are 2x2 matrices: 

f= (^o)' PmN= \7fSI)' DMf= {(DuTy D f)- (R2) 

with F$ N = d M A$ - d N A§ and D M T = d M T + i{A { $ - A$)T = d M T + 2iA M T the usual field 
strength and covariant derivative, where the definition (f4.1|) has been substituted. The STr means 
that one has to symmetrice in Fmn, DmT, T after expanding the square root, and then take the 
trace. 

The expression ( [ti.l| ) is quite involved but we will see that in the particular case we are interested, 
one can deal with it: we will consider quadratic excitations of the gauge fields, while the tachyon 
phase is set to a trivial constant and the tachyon modulus is a non-trivial ^-dependent function (so 
we have to keep all orders in r, d z r). We again take a gauge with A$ = 0. This is enough to compute 
the vector and axial spectrum in the tachyon background. 



15 We adapt it to our present framework, for instance defining the covariant derivative with a different sign and 
disregarding the field. With respect to the main text, we will fix the value of some of the constants that we have 



defined, namely g\ = 2ira' , A = l/(2na'), K = 1, /i 2 = 2. Regarding (|3.8|), this implies R 2 = 3/2. This is inessential 
(the constants can be easily restored) and we have done it for the sake of clarity of the equations. Our convention will 
be that indices M, N running over the five space-time coordinates are contracted with the metric gMN whereas indices 
(i. v running over the Minkowski directions are contracted with the flat metric 77 M „. 
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So let us compute the quadratic expansion in gauge fields. There are terms in A 2 coming from 
the covariant derivatives and terms with F^ 2 . In principle, there could be terms with, schematically, 
ird z rAF coming from a DTDT F product, but these terms would make the action complex and are 
removed by the symmetric trace prescription. In the following, we make the computation in two 
steps: we first compute the A 2 ^ terms and then compute the F 2 terms. 

In order to compute the A 2 terms, we can consider the action: 



Sa? = — STr / d A xdz e 



-f 2 



det(g M N + D M TD N f) 



= — J d 4 xdz a/— det g STr 



- t2 -/det« + D M TD N f) 



(H.3) 



We now compute the determinant. Being inside a STr, the T matrices can be considered as commuting 



objects, so ydet« + D M TD N f) = \f\ + D m TDmT . We have to expand the square root. In order 
to simplify notation, let us define Sj, the coefficients of such expansion: 



v / i+e = E(- 1 ) 

j=0 



3 + 



i (2; -3)11 

j!2i 



3=0 



(H.4) 



Thus, also expanding the exponential: 



C A 2 = -y/-detg^2 



k=0 



Sj STr 

j=0 



T 2k (D M TD M ry 



(H.5) 



The next step is to perform the symmetriced trace, and a major simplification comes out because of 
the particular computation we are doing. Define: 



Ji 



1 

1 



J-2 



-1 

1 



TV. 



such that (use T = T* 

f = tJi D M f = 8 m tJ\ - 2iA M rJ 2 
The order zero term in Am, i-e. the action for the tachyon modulus is just 

C T = -2 V- det g e'^ y/l + d M rd M r, 



(H.6) 



(H.7) 



(H.8) 



the same used in the main text. This means that the discussion of section |3] still holds. 

We now isolate the quadratic term in A M . One out of the j factors of (D M TD M T) 3 has to be 
g^ {—2ir J2) 2 A^A^ while the other j — I factors are Ji9zz(d z T) 2 each. Notice there cannot be crossed 
terms because 8mtA m = in the case we are conisdering. Thus: 



C 



A 2 



9xx 



k=0 



E W T 2k (g; z \d z T) 2 y- l (-4A^T 2 )STr 
3=0 



f2k+2j-2 f2 
J 1 J 2 



(H.9) 
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The j at the beginning of the second line of course comes because the A^A^ term can be chosen 
from any of the j factors in (D m TDmTY . In order to perform the symmetriced trace, notice that, 
in general: 

STrfJ^ 2 ] = ((„ + l)Tr[J 2 ] + nTt&ViJ*]) = ~2n~+l (H ' 10) 

where we have used that Jf is the 2x2 identity matrix and Tr[Ji J 2 J\ J2] = — 7V[Jf] = 2. Substituting: 

00 (_-\\k 00 I 

l# = -Zg^v^Tg E -JT^ £ fe 1 ^) 2 )^ 1 ^ ^,.. _ x (h.h) 

fc=0 ' j=0 ^ 

We should now resum the series. Let us use the identity 16 : 

f U y 3 S J y> = V - f - ^ ia (H.12) 



Using x = —r 2 and y = g z J-(d z r) 2 , we finally find: 



-T 2 a 2 

e ' 



C A2 = -Ag- x W-detg t 2 A u A» / - =da (H.13) 

V 1 + gzz(d z r) 2 a 2 



Let us now compute the F 2 terms. We want to expand the determinant of (ELI) to second order in 
F but to all orders in D Z TD Z T. The determinant reads: 



-det(g MN +F MN +D M TD N T) = gUg zz + D z fD z f) + ^F^g 2 xx (g zz + D z fD z f)+g 3 xx F 2 z (H.14) 
and thus the F 2 contribution to the square root is: 



\gh\l (1 + 9zzDzTD z T)F l j t , u F* lv + \g„gj 1 F% (H.15) 

V 1 + 0£D.TD M T 

Let us start by computing the term with F^F^ U . We have to expand before taking the symmetrized 
trace. Notice that now D Z T = d z rJi up to subleading terms. To shorten notation, we define: 

X = -T 2 , y = g 7z \d z T) 2 (H.16) 

The F^ V F^ U term in the lagrangian density (H.l) then reads: 



1 00 k 

C n» = —A E Ti E ^Tr[Ji k+2j F,Jn (H.17) 



k\ 

k=0 j=0 



16 In order to prove this, notice that Yl'jLo3 s jy : ' = zji+y and consider an auxiliary function 17(a) = 

J2Z tf^E^o 2i+2j-i y ja2i+2i ~ 1 such that d °-9( a ) = a ~ 2 e xa2 2 jl + v a -2 1 ■ Since = and si 1 ) is what we want 

to compute, we arrive at ( H.12| ). 
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It is now easy to compute the symmetriced trace: 



STr[jf +2j F^FH 



2k + 2 j + l V k + 3 + + (* + MWi^J (H.18 



Now, TrfF^F^] = F$ + Fffl which, splitting in vector and axial part and using nota- 
tion of section | gives Tr[F liV F ttv ] = 2V 11V V» V + 2A IIV N 1V . Similarly, Tr[J 1 F fa ,j 1 F' u '] = 2F^ ) F^ = 
2V^V^ - 2A ilv A» v and we find Sl^jf +2j F^F^) = 2V^V^ U + 2i+ 2 2j+l A^A^ . Inserting this in 

(MM-- 



^ 00^,00 / 

c n» = E fr E s ^ + 

fc=o ' i=o x 



2z + 2j + l M 



(H.19) 



The fact that for non-trivial tachyon the symmetric trace produces a coupling between the left and 
right gauge fields was already pointed out in |28| . It results in different kinetic terms for vectors and 
axials. We skip the details of the similar computation leading to F 



2 . 

(J.Z ■ 



£f,2 — —9xx9z 



a/1 + a 2 y 



da (d z A^ 



(H.20) 



By comparing ( |H.19|) , (|H.20|) to ( |4.5|) , we find that the quadratic action for the vector excitation 
is identical regardless the choice between Sen's and Garousi's actions. Nevertheless, the axial part 
changes. From ( |H.19| ), ( |H.20| ) it can be read that, introducing notation of appendix |D| 



■'axial 



-B{z)A llv A^ + A(z)(d z A,) 2 + M{z)A 



(H.21) 



with: 



A(z) = g xx gj J e- T a (y 1 + a?g£(d z T)*) da , 

B{z) = gl z f e- T2a2 ^/l + a 2 g- z \d z Tfda, 
Jo 

^ xx gir 2 / e-^ Ul + a 2 g- z \d z r) 2 ) da 



M(z 



(H.22) 



In order to proceed further, we need estimate the integrals near the IR, where z — > z\ and the tachyon 
diverges, see section 3J.. We will use that in the limit where — x ^> 1 and y ^> 1 with — | ^> 1, it 
happens that: 



2x 



e xa \^/l + ya 2 )- 1 da 



logfe) - 7 + 21og2 
2^ 



(H.23) 
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where 7 is Euler's constant. The first equality is just found by neglecting the 1 inside the square 
root. For the second computation, it is not possible to directly neglect the 1 since the result would 
be divergent, but one can express the integral as /o(\/l + V a?)~ l da + L(e xa2 — + ya 2 )~ 1 da, 

such that the first integral can be done explicitly and in the second one the 1 inside the square root 
can be neglected. 

We are now ready to compute the leading IR behaviour of the Schrodinger potential which will 
determine the behaviour of the highly excited states. We will use that near the IR (z A - z C 1), we 
have: 



K 2 K 1 z A 



z 



z 2 A 5(z A 
9z A C 2 , _ ± 



T = y/—X ~ C (z A — z) 



' 20 



9»(P*f = V « -^(^ ~ *)~™ (H.24) 



From flH.22|) -( [H.24D one can readily check that lirn^^ M(z)/ B(z) = and therefore the term h(u) 



does not contribute in the IR to the Schrodinger potential ( p.7|) . On the other hand, we can obtain 
the relation of the ^-coordinate to the w-coordinate of the Schrodinger problem ( p.5|) : 



§4"lf / ¥r ' ««^W-log(6(z A -^)). (H.25) 
A(z) 20 y/-log{b{z A -z)) 10 

where b is a positive constant that will not be important in the following. We also compute: 

E = (AB)* ~ e (H.26) 

where we have not written multiplicative constants and powers of u which do not affect the leading 
IR behaviour of the Schrodinger potential. Finally, from ( |D.7| ) we find 

25 

V{u) « -^u 2 (H.27) 
vz A 

Since we have a quadratic potential in the IR, the behaviour for asymptotically highly excited axials 



is still Regge-like. Unlike in the main text - section [4.2.1| -, the slope found using Garousi's action 
does not depend on the constant k. Comparing to the vector modes - section |4.1.1| -, we see that the 
Regge slope for the axials is slightly larger, in particular A\/Ay = 10/9, where Ay t A are defined as 
in Q. 
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